(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 11.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     80615,       1981]
NotebookOptionsPosition[     72324,       1863]
NotebookOutlinePosition[     72677,       1879]
CellTagsIndexPosition[     72634,       1876]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[TextData[{
 StyleBox["\nThe program ",
  FontWeight->"Bold",
  FontSlant->"Italic"],
 StyleBox["Motions",
  FontSlant->"Italic"],
 StyleBox["\n",
  FontWeight->"Bold",
  FontSlant->"Italic"]
}], "Subtitle",
 Evaluatable->False,
 CellChangeTimes->{{3.722761415224468*^9, 3.7227614225330725`*^9}, 
   3.722763477638951*^9, 3.7227642159331713`*^9, 3.723978725460906*^9},
 TextAlignment->Center,
 AspectRatioFixed->True,
 FontFamily->"Times New Roman",
 FontSize->36,
 FontWeight->"Bold",ExpressionUUID->"df9343b8-c2ab-43c3-8555-2ada59357c0d"],

Cell[TextData[{
 StyleBox["The programs have been written with ",
  FontVariations->{"CompatibilityType"->0}],
 StyleBox["Mathematica",
  FontSlant->"Italic",
  FontVariations->{"CompatibilityType"->0}],
 StyleBox[" 11. \nOriginal script from A. Romano, Classical Mechanics with \
Mathematica, \nModeling and Simulation in Science, Engineering and \
Technology, \nDOI 10.1007/978-0-8176-8352-8 1,\n\[Copyright] Springer \
Science+Business Media New York 2012\nModified in some part by SeB",
  FontVariations->{"CompatibilityType"->0}]
}], "Text",
 CellChangeTimes->{
  3.487296773359375*^9, {3.487296851765625*^9, 3.487297335328125*^9}, {
   3.488516516809531*^9, 3.488516587606406*^9}, {3.53363855521875*^9, 
   3.533638561953125*^9}, 3.533639119484375*^9, {3.533639440328125*^9, 
   3.533639455171875*^9}, {3.70652439955762*^9, 3.706524408106931*^9}, {
   3.7227562206959224`*^9, 3.7227562212942553`*^9}, 3.7227563039072022`*^9, {
   3.722761513463554*^9, 3.7227615339785776`*^9}, 3.722761622107834*^9, {
   3.723968573667286*^9, 3.7239686154855795`*^9}, {3.7240483797241554`*^9, 
   3.724048412656519*^9}, 3.7240484686770153`*^9, {3.724048513884387*^9, 
   3.724048529416775*^9}, {3.7240485855923634`*^9, 3.724048587221694*^9}, {
   3.7240486213233175`*^9, 3.724048764086769*^9}, 3.7240490037224803`*^9, 
   3.724053537247161*^9, {3.7240535707469463`*^9, 3.72405359365067*^9}, {
   3.724053853279383*^9, 3.7240538697763414`*^9}, {3.724057384625467*^9, 
   3.7240575142497787`*^9}, {3.7240588116837683`*^9, 
   3.7240588549731255`*^9}, {3.724081640472859*^9, 3.724081660675167*^9}, {
   3.724129649026118*^9, 3.7241296770743227`*^9}, {3.7262058150918217`*^9, 
   3.726205839810589*^9}, {3.802103861390992*^9, 3.802103944289927*^9}},
 FontSize->12,ExpressionUUID->"9e441257-d59f-4c27-9fa8-64ea2349ba29"],

Cell[CellGroupData[{

Cell[TextData[{
 StyleBox["Aim of the program ", "Section",
  FontFamily->"Times New Roman",
  FontSize->18],
 StyleBox["Motions",
  FontWeight->"Bold",
  FontSlant->"Italic"]
}], "Section",
 CellChangeTimes->{
  3.7227614629326544`*^9, {3.722761724598725*^9, 3.722761741429019*^9}, 
   3.723978733900426*^9},
 FontFamily->"Times New Roman",
 FontSize->18,ExpressionUUID->"2e9c14ac-be8d-4a96-9d34-cc3f155d6741"],

Cell[TextData[{
 StyleBox["The program ",
  FontVariations->{"CompatibilityType"->0}],
 StyleBox["Motions",
  FontWeight->"Bold",
  FontSlant->"Italic"],
 StyleBox[" ",
  FontVariations->{"CompatibilityType"->0}],
 "analyzes the motion of a material point acted upon by a force F. The \
program provides the parametric equations of motion and the trajectory."
}], "Text",
 CellChangeTimes->{
  3.487296773359375*^9, {3.487296851765625*^9, 3.487297335328125*^9}, {
   3.488516516809531*^9, 3.488516587606406*^9}, {3.53363855521875*^9, 
   3.533638561953125*^9}, 3.533639119484375*^9, {3.533639440328125*^9, 
   3.533639455171875*^9}, {3.70652439955762*^9, 3.706524408106931*^9}, {
   3.7227562206959224`*^9, 3.7227562212942553`*^9}, 3.7227563039072022`*^9, {
   3.722761513463554*^9, 3.7227615339785776`*^9}, 3.722761622107834*^9, {
   3.723978751470697*^9, 3.7239787564218483`*^9}, {3.7239789677884274`*^9, 
   3.7239789892346306`*^9}, {3.723979033533309*^9, 3.723979035834738*^9}},
 FontSize->12,ExpressionUUID->"7c93fce4-aeec-4fcd-8316-27364d5517c4"]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
 StyleBox["Input data of the program ", "Section",
  FontFamily->"Times New Roman",
  FontSize->18],
 StyleBox["Motions",
  FontWeight->"Bold",
  FontSlant->"Italic"]
}], "Section",
 CellChangeTimes->{
  3.7227614629326544`*^9, {3.722761718751762*^9, 3.7227617437109003`*^9}, 
   3.7239787646766577`*^9},
 FontFamily->"Times New Roman",
 FontSize->18,ExpressionUUID->"092d8b89-2f67-40c9-a994-67dc31fd59a8"],

Cell["The command raw of the program is", "Text",
 CellChangeTimes->{
  3.487296773359375*^9, {3.487296851765625*^9, 3.487297335328125*^9}, {
   3.488516516809531*^9, 3.488516587606406*^9}, {3.53363855521875*^9, 
   3.533638561953125*^9}, 3.533639119484375*^9, {3.533639440328125*^9, 
   3.533639455171875*^9}, {3.70652439955762*^9, 3.706524408106931*^9}, {
   3.7227562206959224`*^9, 3.7227562212942553`*^9}, 3.7227563039072022`*^9, {
   3.722761513463554*^9, 3.7227615339785776`*^9}, 3.722761622107834*^9, {
   3.722761758242217*^9, 3.722761764828006*^9}, {3.722761990184806*^9, 
   3.722762075228257*^9}, {3.7233845613088675`*^9, 3.7233845653870535`*^9}},
 FontSize->12,ExpressionUUID->"04780f57-f191-4331-8c7c-9798e27441d0"],

Cell[BoxData[
 FormBox[
  StyleBox[
   RowBox[{
    StyleBox["Motions",
     FontSlant->"Italic"], "[", 
    RowBox[{
    "m", ",", "unk", ",", "F", ",", "P0", ",", "V0", ",", "T1", ",", "T2", 
     ",", "steps", ",", " ", "tf"}], "]"}],
   FontWeight->"Bold"], TraditionalForm]], "DisplayFormula",
 CellChangeTimes->{{3.7227615797960277`*^9, 3.722761651939048*^9}, 
   3.722761692943643*^9, 3.7227620776741457`*^9, {3.7239787683786664`*^9, 
   3.7239787893252954`*^9}},
 TextAlignment->Center,
 FontSize->12,
 FontSlant->"Italic",ExpressionUUID->"844f2f14-e18e-4873-887a-4bcb4ce7f67b"],

Cell["where", "Text",
 CellChangeTimes->{
  3.487296773359375*^9, {3.487296851765625*^9, 3.487297335328125*^9}, {
   3.488516516809531*^9, 3.488516587606406*^9}, {3.53363855521875*^9, 
   3.533638561953125*^9}, 3.533639119484375*^9, {3.533639440328125*^9, 
   3.533639455171875*^9}, {3.70652439955762*^9, 3.706524408106931*^9}, {
   3.7227562206959224`*^9, 3.7227562212942553`*^9}, 3.7227563039072022`*^9, {
   3.722761513463554*^9, 3.7227615339785776`*^9}, 3.722761622107834*^9, {
   3.722761758242217*^9, 3.722761764828006*^9}, {3.722761990184806*^9, 
   3.722762090180395*^9}},
 FontSize->12,ExpressionUUID->"88a7d7c6-3f2f-4a34-8a1b-10caf5a44030"],

Cell[CellGroupData[{

Cell["m = mass of P;", "Item",
 CellChangeTimes->{{3.7227615797960277`*^9, 3.722761651939048*^9}, 
   3.722761692943643*^9, {3.7227617867300777`*^9, 3.722761803082157*^9}, {
   3.7239788015647993`*^9, 3.7239788070669975`*^9}},
 FontSize->12,ExpressionUUID->"96ebc100-b7a2-4c58-9483-a751106aca3e"],

Cell["\<\
unk = list of the coordinate functions (the length of list is 2 or 3 \
according to whether the motion of P takes place in the plane or in the \
space);\
\>", "Item",
 CellChangeTimes->{{3.7227615797960277`*^9, 3.722761651939048*^9}, 
   3.722761692943643*^9, {3.7227617867300777`*^9, 3.722761821634748*^9}, {
   3.723978831213931*^9, 3.7239788377556906`*^9}},
 FontSize->12,ExpressionUUID->"ebd1d692-89e7-4ff8-8091-e01e2f1ce9bd"],

Cell["F = acting force;", "Item",
 CellChangeTimes->{{3.7227615797960277`*^9, 3.722761651939048*^9}, 
   3.722761692943643*^9, {3.7227617867300777`*^9, 3.7227618354129677`*^9}, {
   3.723978846698559*^9, 3.7239788495248833`*^9}},
 FontSize->12,ExpressionUUID->"0428a09f-be09-4e1c-8ffd-0c577684066b"],

Cell["P0 = initial position;", "Item",
 CellChangeTimes->{{3.7227615797960277`*^9, 3.722761651939048*^9}, 
   3.722761692943643*^9, {3.7227617867300777`*^9, 3.7227618383483295`*^9}, {
   3.723978852858414*^9, 3.7239788588530235`*^9}},
 FontSize->12,ExpressionUUID->"800b51e4-5a10-4045-ab95-b437f2cd7467"],

Cell["V0 = initial velocity;", "Item",
 CellChangeTimes->{{3.7227615797960277`*^9, 3.722761651939048*^9}, 
   3.722761692943643*^9, {3.7227617867300777`*^9, 3.7227618418820252`*^9}, {
   3.723978866746561*^9, 3.7239788705148067`*^9}},
 FontSize->12,ExpressionUUID->"48f54f2a-2bbe-4ae4-a715-924236549001"],

Cell["(T1, T2) = time interval in which the motion is analyzed;", "Item",
 CellChangeTimes->{{3.7227615797960277`*^9, 3.722761651939048*^9}, 
   3.722761692943643*^9, {3.7227617867300777`*^9, 3.722761847227234*^9}, 
   3.723978894428373*^9},
 FontSize->12,ExpressionUUID->"36b72aa8-e58d-4a61-a4d9-e4760dfef2a2"],

Cell["steps = steps of numerical integration;", "Item",
 CellChangeTimes->{{3.7227615797960277`*^9, 3.722761651939048*^9}, 
   3.722761692943643*^9, {3.7227617867300777`*^9, 3.722761850866785*^9}},
 FontSize->12,ExpressionUUID->"411034f7-dc86-44ef-8bf4-6c79bdd98508"],

Cell[TextData[{
 StyleBox["tf = ",
  FontColor->GrayLevel[0]],
 "is 0 if ",
 Cell[BoxData[
  FormBox["F", TraditionalForm]],ExpressionUUID->
  "40106d26-16fc-45aa-a305-fde1004f54bc"],
 " is not Newtonian and 1 if it is Newtonian",
 StyleBox[".",
  FontColor->GrayLevel[0]]
}], "Item",
 CellChangeTimes->{{3.7227615797960277`*^9, 3.722761651939048*^9}, 
   3.722761692943643*^9, {3.7227617867300777`*^9, 3.7227618693504725`*^9}, {
   3.723384717686514*^9, 3.723384718656189*^9}, {3.723978908859836*^9, 
   3.7239789392041006`*^9}, {3.80141208065658*^9, 3.801412092914225*^9}},
 FontSize->12,ExpressionUUID->"2823ec98-62bb-4cc1-88a4-c559fbecac7a"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
 StyleBox["The program ", "Section",
  FontFamily->"Times New Roman",
  FontSize->18],
 StyleBox["Motions",
  FontWeight->"Bold",
  FontSlant->"Italic"]
}], "Section",
 CellChangeTimes->{
  3.7227614629326544`*^9, {3.722761718751762*^9, 3.7227617437109003`*^9}, 
   3.722762321053281*^9, 3.723979135301016*^9},
 FontFamily->"Times New Roman",
 FontSize->18,ExpressionUUID->"be9ad61f-e973-403d-84b4-e561c54b767e"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Motions", "[", 
   RowBox[{
   "m_", ",", "unk_", ",", "F_", ",", "P0_", ",", "V0_", ",", "T1_", ",", 
    "T2_", ",", "steps_", ",", " ", "tf_"}], "]"}], ":=", "\n", "\t", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"(*", 
      RowBox[{
      "dim", ",", "U", ",", "en", ",", "sost", ",", "sys", ",", "sol", ",", 
       "m1", ",", "c", ",", "pplot"}], "*)"}], "}"}], ",", "\n", 
    RowBox[{
     RowBox[{"dim", "=", 
      RowBox[{"Length", "[", "unk", "]"}]}], ";", "\n", 
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{"tf", "==", "0"}], ",", 
       RowBox[{"Goto", "[", "2", "]"}], ",", 
       RowBox[{"Goto", "[", "1", "]"}]}], "]"}], ";", "\n", 
     RowBox[{"Label", "[", "1", "]"}], ";", "\n", 
     RowBox[{"U", "=", 
      RowBox[{
       RowBox[{"Integrate", "[", 
        RowBox[{"F", ",", "x"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}], ";", "\n", 
     RowBox[{"en", "=", 
      RowBox[{
       RowBox[{
        RowBox[{
         FractionBox["1", "2"], "m", " ", 
         RowBox[{"Sum", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"(", 
             RowBox[{"V0", "[", 
              RowBox[{"[", "i", "]"}], "]"}], ")"}], "^", "2"}], ",", 
           RowBox[{"{", 
            RowBox[{"i", ",", "1", ",", "dim"}], "}"}]}], "]"}]}], "-", "U"}],
        "/.", 
       RowBox[{"Table", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"unk", "[", 
           RowBox[{"[", "i", "]"}], "]"}], "\[Rule]", 
          RowBox[{"P0", "[", 
           RowBox[{"[", "i", "]"}], "]"}]}], ",", 
         RowBox[{"{", 
          RowBox[{"i", ",", "1", ",", "dim"}], "}"}]}], "]"}]}]}], ";", "\n", 
     RowBox[{"Print", "[", 
      RowBox[{
       RowBox[{"Style", "[", 
        RowBox[{"\"\<Total energy = \>\"", ",", "\"\<Output\>\"", ",", 
         RowBox[{"FontFamily", "\[Rule]", "\"\<Times-Plain\>\""}]}], "]"}], 
       ",", 
       RowBox[{"Style", "[", 
        RowBox[{"en", ",", "\"\<Output\>\"", ",", 
         RowBox[{"FontFamily", "\[Rule]", "\"\<Times-Plain\>\""}]}], "]"}]}], 
      "]"}], ";", "\n", 
     RowBox[{"Label", "[", "2", "]"}], ";", "\[IndentingNewLine]", "\n", 
     RowBox[{"(*", " ", 
      RowBox[{"Equations", " ", "of", " ", "motion"}], " ", "*)"}], "\n", 
     "\t\t", 
     RowBox[{"sost", "=", 
      RowBox[{
       RowBox[{"Table", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{
            RowBox[{"unk", "[", 
             RowBox[{"[", "i", "]"}], "]"}], "->", 
            RowBox[{
             RowBox[{"unk", "[", 
              RowBox[{"[", "i", "]"}], "]"}], "[", "t", "]"}]}], ",", 
           RowBox[{
            RowBox[{
             RowBox[{"unk", "[", 
              RowBox[{"[", "i", "]"}], "]"}], "'"}], "->", 
            RowBox[{
             RowBox[{
              RowBox[{"unk", "[", 
               RowBox[{"[", "i", "]"}], "]"}], "'"}], "[", "t", "]"}]}], ",", 
           "\n", "\t\t\t\t\t", 
           RowBox[{
            RowBox[{
             RowBox[{
              RowBox[{"unk", "[", 
               RowBox[{"[", "i", "]"}], "]"}], "'"}], "'"}], "->", 
            RowBox[{
             RowBox[{
              RowBox[{
               RowBox[{"unk", "[", 
                RowBox[{"[", "i", "]"}], "]"}], "'"}], "'"}], "[", "t", 
             "]"}]}]}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{"i", ",", "1", ",", "dim"}], "}"}]}], "]"}], "//", 
       "Flatten"}]}], ";", "\n", "\t\t", 
     RowBox[{"sys1", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"m", " ", 
          RowBox[{
           RowBox[{
            RowBox[{"unk", "[", 
             RowBox[{"[", "i", "]"}], "]"}], "'"}], "'"}]}], "==", 
         RowBox[{"F", "[", 
          RowBox[{"[", "i", "]"}], "]"}]}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "dim"}], "}"}]}], "]"}]}], ";", "\n", 
     "\t\t", 
     RowBox[{"sys", "=", 
      RowBox[{"sys1", "/.", "sost"}]}], ";", "\n", 
     RowBox[{"(*", 
      RowBox[{"Initial", " ", "data"}], "*)"}], "\n", "\t\t", 
     RowBox[{"indata", "=", 
      RowBox[{
       RowBox[{"Table", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{
            RowBox[{
             RowBox[{"unk", "[", 
              RowBox[{"[", "j", "]"}], "]"}], "[", "0", "]"}], "\[Equal]", 
            RowBox[{"P0", "[", 
             RowBox[{"[", "j", "]"}], "]"}]}], ",", 
           RowBox[{
            RowBox[{
             RowBox[{
              RowBox[{"unk", "[", 
               RowBox[{"[", "j", "]"}], "]"}], "'"}], "[", "0", "]"}], 
            "\[Equal]", 
            RowBox[{"V0", "[", 
             RowBox[{"[", "j", "]"}], "]"}]}]}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{"j", ",", "1", ",", "dim"}], "}"}]}], "]"}], "//", 
       "Flatten"}]}], ";", "\n", "\[IndentingNewLine]", 
     RowBox[{"(*", 
      RowBox[{"Numeric", " ", "solution"}], "*)"}], "\n", "\t", 
     RowBox[{"sol", "=", 
      RowBox[{"NDSolve", "[", 
       RowBox[{
        RowBox[{"Join", "[", 
         RowBox[{"sys", ",", "indata"}], "]"}], ",", "unk", ",", 
        RowBox[{"{", 
         RowBox[{"t", ",", "T1", ",", "T2"}], "}"}], ",", 
        RowBox[{"MaxSteps", "->", "steps"}]}], "]"}]}], ";", "\n", 
     "\[IndentingNewLine]", 
     RowBox[{"(*", 
      RowBox[{"Coordinates", " ", "versus", " ", "time"}], "*)"}], "\n", 
     RowBox[{"m1", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"unk", "[", 
         RowBox[{"[", "i", "]"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "dim"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"color", "=", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"RGBColor", "[", 
         RowBox[{"0", ",", "0", ",", "0"}], "]"}], ",", 
        RowBox[{"RGBColor", "[", 
         RowBox[{"1", ",", "0", ",", "0"}], "]"}], ",", 
        RowBox[{"RGBColor", "[", 
         RowBox[{"0", ",", "0", ",", "1"}], "]"}]}], "}"}]}], ";", 
     "\[IndentingNewLine]", "\t\t", 
     RowBox[{
      RowBox[{"c", "[", "i_", "]"}], ":=", 
      RowBox[{"Plot", "[", 
       RowBox[{
        RowBox[{"Evaluate", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"unk", "[", 
            RowBox[{"[", "i", "]"}], "]"}], "[", "t", "]"}], "/.", "sol"}], 
         "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"t", ",", "T1", ",", "T2"}], "}"}], ",", 
        RowBox[{"AxesLabel", "->", 
         RowBox[{"{", 
          RowBox[{"t", ",", "\"\<\>\""}], "}"}]}], ",", 
        RowBox[{"PlotLegends", "\[Rule]", 
         RowBox[{"LineLegend", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"color", "[", 
             RowBox[{"[", "i", "]"}], "]"}], "}"}], ",", 
           RowBox[{"{", 
            RowBox[{"unk", "[", 
             RowBox[{"[", "i", "]"}], "]"}], "}"}]}], "]"}]}], ",", 
        RowBox[{"PlotRange", "->", "All"}], ",", 
        RowBox[{"PlotStyle", "\[Rule]", 
         RowBox[{"color", "[", 
          RowBox[{"[", "i", "]"}], "]"}]}]}], "]"}]}], ";", "\n", "\t\t", 
     RowBox[{"pplot", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"c", "[", "i", "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "dim"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"StylePrint", "[", 
      RowBox[{"\"\<Equation of motion\>\"", ",", "\"\<Output\>\"", ",", 
       RowBox[{"FontFamily", "\[Rule]", "\"\<Times-Plain\>\""}]}], "]"}], ";",
      "\n", "\t\t", 
     RowBox[{"Print", "[", 
      RowBox[{"Show", "[", "pplot", "]"}], "]"}], ";", "\n", "\n", 
     RowBox[{"(*", 
      RowBox[{
       RowBox[{
        RowBox[{"Trajectories", " ", "for", " ", "n"}], "=", "2"}], ",", 
       "3"}], "*)"}], "\[IndentingNewLine]", "\t\t", 
     RowBox[{"Which", "[", 
      RowBox[{
       RowBox[{"dim", "==", "1"}], ",", 
       RowBox[{"Goto", "[", "5", "]"}], ",", 
       RowBox[{"dim", "==", "2"}], ",", 
       RowBox[{"Goto", "[", "3", "]"}], ",", 
       RowBox[{"dim", "==", "3"}], ",", 
       RowBox[{"Goto", "[", "4", "]"}]}], "]"}], ";", "\n", "\t\t", 
     RowBox[{"Label", "[", "3", "]"}], ";", "\[IndentingNewLine]", 
     RowBox[{"StylePrint", "[", 
      RowBox[{"\"\<Trajectory\>\"", ",", "\"\<Output\>\"", ",", 
       RowBox[{"FontFamily", "\[Rule]", "\"\<Times-Plain\>\""}]}], "]"}], ";",
      "\n", "\t\t", 
     RowBox[{"Print", "[", 
      RowBox[{"ParametricPlot", "[", 
       RowBox[{
        RowBox[{"Evaluate", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{
             RowBox[{"unk", "[", 
              RowBox[{"[", "1", "]"}], "]"}], "[", "t", "]"}], ",", 
            RowBox[{
             RowBox[{"unk", "[", 
              RowBox[{"[", "2", "]"}], "]"}], "[", "t", "]"}]}], "}"}], "/.", 
          "sol"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"t", ",", "T1", ",", "T2"}], "}"}], ",", 
        RowBox[{"AxesLabel", "->", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"unk", "[", 
            RowBox[{"[", "1", "]"}], "]"}], ",", 
           RowBox[{"unk", "[", 
            RowBox[{"[", "2", "]"}], "]"}]}], "}"}]}], ",", 
        RowBox[{"PlotRange", "->", "All"}], ",", 
        RowBox[{"PlotStyle", "\[Rule]", "Black"}]}], "]"}], "]"}], ";", "\n", 
     "\t\t", 
     RowBox[{"Goto", "[", "5", "]"}], ";", "\n", "\t", 
     RowBox[{"Label", "[", "4", "]"}], ";", "\[IndentingNewLine]", 
     RowBox[{"StylePrint", "[", 
      RowBox[{"\"\<Trajectory\>\"", ",", "\"\<Output\>\"", ",", 
       RowBox[{"FontFamily", "\[Rule]", "\"\<Times-Plain\>\""}]}], "]"}], ";",
      "\n", "\t\t\t\t", 
     RowBox[{"Print", "[", 
      RowBox[{"ParametricPlot3D", "[", 
       RowBox[{
        RowBox[{"Evaluate", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{
             RowBox[{"unk", "[", 
              RowBox[{"[", "1", "]"}], "]"}], "[", "t", "]"}], ",", 
            RowBox[{
             RowBox[{"unk", "[", 
              RowBox[{"[", "2", "]"}], "]"}], "[", "t", "]"}], ",", 
            RowBox[{
             RowBox[{"unk", "[", 
              RowBox[{"[", "3", "]"}], "]"}], "[", "t", "]"}]}], "}"}], "/.", 
          "sol"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"t", ",", "T1", ",", "T2"}], "}"}], ",", 
        RowBox[{"AxesLabel", "->", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"unk", "[", 
            RowBox[{"[", "1", "]"}], "]"}], ",", 
           RowBox[{"unk", "[", 
            RowBox[{"[", "2", "]"}], "]"}], ",", 
           RowBox[{"unk", "[", 
            RowBox[{"[", "3", "]"}], "]"}]}], "}"}]}], ",", 
        RowBox[{"PlotPoints", "->", "150"}], ",", 
        RowBox[{"PlotStyle", "\[Rule]", "Black"}]}], "]"}], "]"}], ";", "\n", 
     "\t", 
     RowBox[{"Label", "[", "5", "]"}], ";"}]}], "\n", "\n", "\t\t", 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.48699314196875*^9, 3.486993151640625*^9}, {
   3.487226727758194*^9, 3.487226756914444*^9}, {3.487297420296875*^9, 
   3.487297490515625*^9}, {3.48826032959375*^9, 3.488260337140625*^9}, {
   3.488516508356406*^9, 3.488516513575156*^9}, {3.488516982825156*^9, 
   3.488517033106406*^9}, 3.488517114684531*^9, {3.488517483590781*^9, 
   3.488517490497031*^9}, {3.488518767059531*^9, 3.488518771012656*^9}, {
   3.488519179403281*^9, 3.488519184043906*^9}, {3.488519233090781*^9, 
   3.488519288981406*^9}, {3.488680561765625*^9, 3.48868056540625*^9}, {
   3.548557899640625*^9, 3.5485579043125*^9}, 3.7227558620892477`*^9, {
   3.7227558933151007`*^9, 3.72275589610955*^9}, 3.722762640540414*^9, 
   3.722762686969385*^9, {3.72276350540983*^9, 3.7227635475629787`*^9}, 
   3.723384751275381*^9, {3.723384788775066*^9, 3.7233848202099047`*^9}, {
   3.7233848955385294`*^9, 3.7233849862093315`*^9}, {3.7233850433699036`*^9, 
   3.7233850489269123`*^9}, 3.7233851696267805`*^9, {3.7233852913764286`*^9, 
   3.7233852979783196`*^9}, {3.723385346969463*^9, 3.723385372063978*^9}, {
   3.7233854235869064`*^9, 3.723385424140579*^9}, {3.7233854819065886`*^9, 
   3.723385482527684*^9}, {3.7233855508282022`*^9, 3.723385607519827*^9}, {
   3.723385727024231*^9, 3.7233857315722685`*^9}, 3.7239791275065804`*^9, 
   3.7239791844514823`*^9, {3.7239795074644985`*^9, 3.723979521256486*^9}, {
   3.7239795613819504`*^9, 3.723979641450824*^9}, {3.7239796734191895`*^9, 
   3.723979673615765*^9}, {3.7239813141308856`*^9, 3.723981320728196*^9}, {
   3.723981401550949*^9, 3.723981415821499*^9}, {3.723982086904264*^9, 
   3.723982090125343*^9}, {3.723982389826375*^9, 3.7239824667885065`*^9}, {
   3.7239825560472097`*^9, 3.7239825564208*^9}, 3.7239826472126265`*^9, {
   3.7239827854088545`*^9, 3.723982795723876*^9}, 3.7239829124733067`*^9, {
   3.7239832454159756`*^9, 3.7239832702436867`*^9}, {3.7239833057913322`*^9, 
   3.7239833380830097`*^9}},
 FontSize->12,
 Background->GrayLevel[0.843763],
 CellLabel->"In[1]:=",ExpressionUUID->"626c1585-ac06-43dc-bf9d-f5a751a84769"]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[StyleBox["Worked Examples",
 FontFamily->"Times New Roman"]], "Section",
 FontSize->18,ExpressionUUID->"9d1eb348-1a14-4b6a-853a-7747f90d2407"],

Cell[CellGroupData[{

Cell[TextData[{
 "E",
 StyleBox["xample",
  FontFamily->"Times New Roman"],
 " seb1: ",
 StyleBox["Material point moving along an inclined plain on straight line and \
acted upon by g1*sin(\[Alpha]).",
  FontWeight->"Bold"]
}], "Subsection",
 Evaluatable->False,
 CellChangeTimes->{{3.7227623633962064`*^9, 3.7227623643477526`*^9}, 
   3.7227624245276814`*^9, 3.723979406813022*^9, 3.7258609339668283`*^9, {
   3.801412158862597*^9, 3.801412232404232*^9}, 3.802149505756115*^9, 
   3.802450346766377*^9},
 AspectRatioFixed->True,
 FontFamily->"Times New Roman",
 FontSize->16,ExpressionUUID->"e7d99775-53c8-443d-bc9b-999bb355f654"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"m1", "=", "1"}], ";"}], "\n", 
 RowBox[{
  RowBox[{"g1", "=", " ", "9.8"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Alpha]1", "=", 
   RowBox[{"\[Pi]", "/", "6"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"unk", "=", 
    RowBox[{"{", "x", "}"}]}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{
   "here", " ", "x", " ", "is", " ", "assumed", " ", "along", " ", "the", " ",
     "slope"}], " ", "*)"}]}], "\n", 
 RowBox[{
  RowBox[{"F", "=", 
   RowBox[{"{", 
    RowBox[{"m1", "*", "g1", "*", 
     RowBox[{"Sin", "[", "\[Alpha]1", "]"}]}], "}"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"P0", "=", 
    RowBox[{"{", "0", "}"}]}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{"initial", " ", "position"}], "*)"}]}], "\n", 
 RowBox[{
  RowBox[{"V0", "=", 
   RowBox[{"{", "0", "}"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"T1", "=", "0"}], ";"}], "\n", 
 RowBox[{
  RowBox[{"T2", "=", "10"}], ";"}], "\n", 
 RowBox[{
  RowBox[{"steps", "=", "5000"}], ";"}], "\n", 
 RowBox[{
  RowBox[{"tf", "=", "0"}], ";"}], "\n", 
 RowBox[{"Motions", "[", 
  RowBox[{
  "m1", ",", "unk", ",", "F", ",", "P0", ",", "V0", ",", "T1", ",", "T2", ",",
    "steps", ",", "tf"}], "]"}]}], "Input",
 CellChangeTimes->{{3.488516597247031*^9, 3.488516600497031*^9}, 
   3.533638394234375*^9, 3.53363920596875*^9, 3.533639427546875*^9, 
   3.533639468875*^9, {3.7227626044903383`*^9, 3.72276261200213*^9}, {
   3.723384253628869*^9, 3.723384294709807*^9}, 3.7239794232175193`*^9, {
   3.80141233142813*^9, 3.801412456743663*^9}, {3.801412945428204*^9, 
   3.801412946968837*^9}, {3.80141394002007*^9, 3.801413961054491*^9}, {
   3.802450198361515*^9, 3.8024502215201073`*^9}},
 FontSize->12,
 CellLabel->"In[2]:=",ExpressionUUID->"bbb28080-82de-4f03-b226-4247ac15aab6"]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
 "E",
 StyleBox["xample",
  FontFamily->"Times New Roman"],
 " seb2: ",
 StyleBox["Material point moving along an inclined plain on straight line and \
acted upon by g1*sin(\[Alpha]) \:0438 \:0440\:0435\:0430\:043a\:0446\:0438\
\:0435\:0439 \:043e\:043f\:043e\:0440\:044b, \:0442\:0440\:0435\:043d\:0438\
\:044f \:043d\:0435\:0442.",
  FontWeight->"Bold"]
}], "Subsection",
 Evaluatable->False,
 CellChangeTimes->{{3.7227623633962064`*^9, 3.7227623643477526`*^9}, {
   3.7227624078391485`*^9, 3.722762412550025*^9}, 3.7239814420362844`*^9, {
   3.8014125894411163`*^9, 3.80141259052518*^9}, {3.801412731179865*^9, 
   3.801412742256505*^9}, {3.801499374901988*^9, 3.801499419838655*^9}, {
   3.802149380097231*^9, 3.802149382699342*^9}, 3.802149510487389*^9, 
   3.802450350123867*^9},
 AspectRatioFixed->True,
 FontFamily->"Times New Roman",
 FontSize->16,ExpressionUUID->"e57fe7c6-2316-45d0-9a5d-1efa01ac6413"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"m1", "=", "1."}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"g1", "=", " ", "9.8"}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{"take", " ", "approximate", " ", "for", " ", "simplicity"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Alpha]1", "=", 
   RowBox[{"\[Pi]", "/", "6"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"unk", "=", 
    RowBox[{"{", 
     RowBox[{"x", ",", "y"}], "}"}]}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"here", " ", "x", " ", "and", " ", "y"}], " ", "-", " ", 
    RowBox[{"usual", " ", "Cartesian"}]}], " ", "*)"}]}], "\n", 
 RowBox[{
  RowBox[{"F", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"m1", "*", "g1", "*", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"Cos", "[", "\[Alpha]1", "]"}], "*", 
        RowBox[{"Sin", "[", "\[Alpha]1", "]"}]}], ")"}]}], ",", 
     RowBox[{"m1", "*", "g1", "*", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{"Cos", "[", "\[Alpha]1", "]"}], ")"}], "^", "2"}], "-", 
        "1"}], ")"}]}]}], "}"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"t1", "=", "0"}], ";"}], "\n", 
 RowBox[{
  RowBox[{"t2", "=", "3"}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"P0", "=", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"0", ",", 
       RowBox[{
        RowBox[{"-", "g1"}], "*", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{"Cos", "[", "\[Alpha]1", "]"}], ")"}], "^", "2"}], "-", 
          "1"}], ")"}]}]}], "}"}], "*", 
     RowBox[{
      RowBox[{"t2", "^", "2"}], "/", "2"}]}]}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{
   "adjusting", " ", "altitude", " ", "to", " ", "falling", " ", "to", " ", 
    "zero"}], " ", "*)"}]}], "\n", 
 RowBox[{
  RowBox[{"V0", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0"}], "}"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"steps", "=", "50"}], ";"}], "\n", 
 RowBox[{
  RowBox[{"tf", "=", "0"}], ";"}], "\n", 
 RowBox[{"Motions", "[", 
  RowBox[{
  "m1", ",", "unk", ",", "F", ",", "P0", ",", "V0", ",", "t1", ",", "t2", ",",
    "steps", ",", "tf"}], "]"}]}], "Input",
 CellChangeTimes->{{3.488516603965781*^9, 3.488516609450156*^9}, 
   3.533638409953125*^9, {3.533639196609375*^9, 3.533639200109375*^9}, 
   3.533639232609375*^9, 3.7227627607774453`*^9, {3.7239814547375064`*^9, 
   3.723981456019337*^9}, {3.8014127751361628`*^9, 3.8014128097557993`*^9}, {
   3.801412840033223*^9, 3.801412840844523*^9}, {3.8014129630674887`*^9, 
   3.801412974090295*^9}, 3.801413298214366*^9, {3.801413365263443*^9, 
   3.8014133711072283`*^9}, {3.801413500442295*^9, 3.801413521750505*^9}, {
   3.801413766388073*^9, 3.8014138548879213`*^9}, {3.80141397735352*^9, 
   3.801413996378757*^9}, {3.8014983672015257`*^9, 3.801498394810019*^9}, {
   3.801498484419816*^9, 3.801498516950513*^9}, {3.801498550119773*^9, 
   3.8014985530588493`*^9}, {3.801498639396269*^9, 3.8014986529798727`*^9}, {
   3.801498684739962*^9, 3.8014987056069527`*^9}, {3.8014987616399307`*^9, 
   3.801498840921955*^9}, 3.801498880172531*^9, {3.801498920206132*^9, 
   3.801498988748966*^9}, {3.801499084421528*^9, 3.801499084634947*^9}, {
   3.801668329589964*^9, 3.801668334450191*^9}, {3.802450241287763*^9, 
   3.8024502570878267`*^9}, {3.802450294646887*^9, 3.802450320876512*^9}},
 FontSize->12,
 CellLabel->"In[14]:=",ExpressionUUID->"157026c9-04c9-4d7a-9418-000cd75bd159"],

Cell[BoxData[" "], "Input",
 CellChangeTimes->{3.801499192671075*^9},
 CellLabel->"In[26]:=",ExpressionUUID->"dc939f53-588f-4af3-bd3c-17fb1dcbb632"]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
 "E",
 StyleBox["xample",
  FontFamily->"Times New Roman"],
 " seb3: ",
 StyleBox["Material point with mass m1 moving along an inclined (angle \
\[Alpha]1) plain (wedge) on straight line and acted upon by gravity and \
reaction of surface. The wedge of mass M1 is moving  on table with \
acceleration a1. Find the motion of m1.",
  FontWeight->"Bold"]
}], "Subsection",
 Evaluatable->False,
 CellChangeTimes->{{3.7227623633962064`*^9, 3.7227623643477526`*^9}, {
   3.7227624078391485`*^9, 3.722762412550025*^9}, 3.7239814420362844`*^9, {
   3.8014125894411163`*^9, 3.80141259052518*^9}, {3.801412731179865*^9, 
   3.801412742256505*^9}, 3.80149934928935*^9, {3.8014994372951393`*^9, 
   3.801499559283052*^9}, {3.801499610578313*^9, 3.801499629038383*^9}, {
   3.801499940257044*^9, 3.8014999526139793`*^9}, {3.801666617251504*^9, 
   3.801666618525972*^9}, {3.8021495159595203`*^9, 3.8021495180469513`*^9}, {
   3.8021498909601717`*^9, 3.802149893126061*^9}, {3.802150056368335*^9, 
   3.80215009336775*^9}, {3.802450358932087*^9, 3.802450382445907*^9}, 
   3.802450895070538*^9},
 AspectRatioFixed->True,
 FontFamily->"Times New Roman",
 FontSize->16,ExpressionUUID->"391a5db5-3c8f-43f7-a11a-6166c9167d30"],

Cell[TextData[{
 "In the first piece we simply check the motion of mass ",
 StyleBox["m1 ",
  FontSlant->"Italic"],
 "on a fixed, non-moving wedge.\nIn this case the reaction is ",
 StyleBox["N1=m1*g1*",
  FontSlant->"Italic"],
 "cos\[Alpha]. Its projections to X and  to Y (together with ",
 StyleBox["m1*g1",
  FontSlant->"Italic"],
 ") gives the force ",
 StyleBox["F1 ",
  FontSlant->"Italic"],
 "below. "
}], "Text",
 CellChangeTimes->{{3.8023337654903316`*^9, 3.802333773134351*^9}, {
  3.8023340700921383`*^9, 3.80233413450178*^9}, {3.802334227816073*^9, 
  3.802334343185988*^9}, {3.8024489409175053`*^9, 3.802448959095964*^9}, {
  3.802450938749468*^9, 
  3.80245097755201*^9}},ExpressionUUID->"8332d20d-e726-4921-9ed8-\
ebe81b7dc649"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"m1", "=", "1."}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"g1", "=", " ", "9.8"}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{"take", " ", "approximate", " ", "for", " ", "simplicity"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Alpha]1", "=", 
   RowBox[{"\[Pi]", "/", "6"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"unk", "=", 
    RowBox[{"{", 
     RowBox[{"x", ",", "y"}], "}"}]}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"here", " ", "x", " ", "and", " ", "y"}], " ", "-", " ", 
    RowBox[{"usual", " ", "Cartesian"}]}], " ", "*)"}]}], "\n", 
 RowBox[{
  RowBox[{"F1", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"m1", "*", "g1", "*", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"Cos", "[", "\[Alpha]1", "]"}], "*", 
        RowBox[{"Sin", "[", "\[Alpha]1", "]"}]}], ")"}]}], ",", 
     RowBox[{"m1", "*", "g1", "*", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{"Cos", "[", "\[Alpha]1", "]"}], ")"}], "^", "2"}], "-", 
        "1"}], ")"}]}]}], "}"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"t1", "=", "0"}], ";"}], "\n", 
 RowBox[{
  RowBox[{"t2", "=", "3"}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"P01", "=", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"0", ",", 
       RowBox[{
        RowBox[{"-", "g1"}], "*", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{"Cos", "[", "\[Alpha]1", "]"}], ")"}], "^", "2"}], "-", 
          "1"}], ")"}]}]}], "}"}], "*", 
     RowBox[{
      RowBox[{"t2", "^", "2"}], "/", "2"}]}]}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{
   "adjusting", " ", "altitude", " ", "to", " ", "falling", " ", "to", " ", 
    "zero"}], " ", "*)"}]}], "\n", 
 RowBox[{
  RowBox[{"V01", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0"}], "}"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"steps", "=", "5000"}], ";"}], "\n", 
 RowBox[{
  RowBox[{"tf", "=", "0"}], ";"}], "\n", 
 RowBox[{"Motions", "[", 
  RowBox[{
  "m1", ",", "unk", ",", "F1", ",", "P0", ",", "V0", ",", "t1", ",", "t2", 
   ",", "steps", ",", "tf"}], "]"}]}], "Input",
 CellChangeTimes->{{3.488516603965781*^9, 3.488516609450156*^9}, 
   3.533638409953125*^9, {3.533639196609375*^9, 3.533639200109375*^9}, 
   3.533639232609375*^9, 3.7227627607774453`*^9, {3.7239814547375064`*^9, 
   3.723981456019337*^9}, {3.8014127751361628`*^9, 3.8014128097557993`*^9}, {
   3.801412840033223*^9, 3.801412840844523*^9}, {3.8014129630674887`*^9, 
   3.801412974090295*^9}, 3.801413298214366*^9, {3.801413365263443*^9, 
   3.8014133711072283`*^9}, {3.801413500442295*^9, 3.801413521750505*^9}, {
   3.801413766388073*^9, 3.8014138548879213`*^9}, {3.80141397735352*^9, 
   3.801413996378757*^9}, {3.8014983672015257`*^9, 3.801498394810019*^9}, {
   3.801498484419816*^9, 3.801498516950513*^9}, {3.801498550119773*^9, 
   3.8014985530588493`*^9}, {3.801498639396269*^9, 3.8014986529798727`*^9}, {
   3.801498684739962*^9, 3.8014987056069527`*^9}, {3.8014987616399307`*^9, 
   3.801498840921955*^9}, 3.801498880172531*^9, {3.801498920206132*^9, 
   3.801498988748966*^9}, {3.801499084421528*^9, 3.801499084634947*^9}, {
   3.801499576164363*^9, 3.8014995763877296`*^9}, {3.801668355011529*^9, 
   3.801668357615589*^9}, {3.8017104340360937`*^9, 3.8017104384798613`*^9}, {
   3.8024503909524612`*^9, 3.8024504439953337`*^9}, {3.802450981961825*^9, 
   3.802451003620372*^9}},
 FontSize->12,
 CellLabel->"In[27]:=",ExpressionUUID->"c68b342b-e3ef-4cb7-b2ab-44b43d21661b"],

Cell["Now consider just the motion of M1 in 1D:", "Text",
 CellChangeTimes->{{3.8024510968925858`*^9, 
  3.802451127623684*^9}},ExpressionUUID->"0b8915e2-a79f-4ff6-b058-\
57f0ca2c773f"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"M1", "=", "3."}], ";"}], "\n", 
 RowBox[{
  RowBox[{"unkM", "=", 
   RowBox[{"{", "x", "}"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"a1", "=", "3."}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{"minus", " ", "if", " ", "moves", " ", "to", " ", "left"}], " ", 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"FM1", "=", 
   RowBox[{"{", 
    RowBox[{"M1", "*", "a1"}], "}"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"tM1", "=", "0"}], ";"}], "\n", 
 RowBox[{
  RowBox[{"tM2", "=", "3"}], ";"}], "\n", 
 RowBox[{
  RowBox[{"PM0", "=", 
   RowBox[{"{", "0", "}"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"VM0", "=", 
   RowBox[{"{", "0", "}"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"steps", "=", "5000"}], ";"}], "\n", 
 RowBox[{
  RowBox[{"tf", "=", "0"}], ";"}], "\n", 
 RowBox[{"Motions", "[", 
  RowBox[{
  "M1", ",", "unkM", ",", "FM1", ",", "PM0", ",", "VM0", ",", "tM1", ",", 
   "tM2", ",", "steps", ",", "tf"}], "]"}]}], "Input",
 CellChangeTimes->{{3.488516603965781*^9, 3.488516609450156*^9}, 
   3.533638409953125*^9, {3.533639196609375*^9, 3.533639200109375*^9}, 
   3.533639232609375*^9, 3.7227627607774453`*^9, {3.7239814547375064`*^9, 
   3.723981456019337*^9}, {3.8014127751361628`*^9, 3.8014128097557993`*^9}, {
   3.801412840033223*^9, 3.801412840844523*^9}, {3.8014129630674887`*^9, 
   3.801412974090295*^9}, 3.801413298214366*^9, {3.801413365263443*^9, 
   3.8014133711072283`*^9}, {3.801413500442295*^9, 3.801413521750505*^9}, {
   3.801413766388073*^9, 3.8014138548879213`*^9}, {3.80141397735352*^9, 
   3.801413996378757*^9}, {3.8014983672015257`*^9, 3.801498394810019*^9}, {
   3.801498484419816*^9, 3.801498516950513*^9}, {3.801498550119773*^9, 
   3.8014985530588493`*^9}, {3.801498639396269*^9, 3.8014986529798727`*^9}, {
   3.801498684739962*^9, 3.8014987056069527`*^9}, {3.8014987616399307`*^9, 
   3.801498840921955*^9}, 3.801498880172531*^9, {3.801498920206132*^9, 
   3.801498988748966*^9}, {3.801499084421528*^9, 3.801499084634947*^9}, {
   3.801499576164363*^9, 3.8014995763877296`*^9}, {3.8014996557479963`*^9, 
   3.801499659715501*^9}, {3.8014996981436043`*^9, 3.8014998754578876`*^9}, {
   3.80166665965838*^9, 3.8016666619636297`*^9}, {3.801666718494713*^9, 
   3.801666742795319*^9}, {3.801668372261248*^9, 3.801668383923596*^9}, {
   3.802450468217289*^9, 3.802450488021326*^9}, 3.802450795088221*^9},
 FontSize->12,
 CellLabel->"In[39]:=",ExpressionUUID->"26e627c5-dd72-4926-9caa-30c7417a47a4"],

Cell["\<\
Consider ref. frame where wedge is at rest. Then it is non-inertial and we \
have equiv. gravity \
\>", "Text",
 CellChangeTimes->{{3.801500034245212*^9, 
  3.801500093033792*^9}},ExpressionUUID->"60d2d6fc-8561-4cd8-886b-\
50aad0651079"],

Cell["\<\
Find acceleration a for M such that m will be at rest. Approximately it is 6, \
analytically g1 tg \[Alpha]1:\
\>", "Text",
 CellChangeTimes->{{3.801674944265684*^9, 3.801674994361802*^9}, {
  3.802104331604205*^9, 3.802104369156797*^9}, {3.8021501087785482`*^9, 
  3.802150114254904*^9}, {3.8024505265177517`*^9, 
  3.8024505282270107`*^9}},ExpressionUUID->"607be239-0850-4601-b0da-\
2ef9ae900605"],

Cell[BoxData[
 RowBox[{"\[Alpha]1", "=", 
  RowBox[{"\[Pi]", "/", "6"}]}]], "Input",
 CellChangeTimes->{{3.488516603965781*^9, 3.488516609450156*^9}, 
   3.533638409953125*^9, {3.533639196609375*^9, 3.533639200109375*^9}, 
   3.533639232609375*^9, 3.7227627607774453`*^9, {3.7239814547375064`*^9, 
   3.723981456019337*^9}, {3.8014127751361628`*^9, 3.8014128097557993`*^9}, {
   3.801412840033223*^9, 3.801412840844523*^9}, {3.8014129630674887`*^9, 
   3.801412974090295*^9}, 3.801413298214366*^9, {3.801413365263443*^9, 
   3.8014133711072283`*^9}, {3.801413500442295*^9, 3.801413521750505*^9}, {
   3.801413766388073*^9, 3.8014138548879213`*^9}, {3.80141397735352*^9, 
   3.801413996378757*^9}, {3.8014983672015257`*^9, 3.801498394810019*^9}, {
   3.801498484419816*^9, 3.801498516950513*^9}, {3.801498550119773*^9, 
   3.8014985530588493`*^9}, {3.801498639396269*^9, 3.8014986529798727`*^9}, {
   3.801498684739962*^9, 3.8014987056069527`*^9}, {3.8014987616399307`*^9, 
   3.801498840921955*^9}, 3.801498880172531*^9, {3.801498920206132*^9, 
   3.801498988748966*^9}, {3.801499084421528*^9, 3.801499084634947*^9}, {
   3.801499576164363*^9, 3.8014995763877296`*^9}, {3.801500116322962*^9, 
   3.8015001235638103`*^9}, {3.801500297090946*^9, 3.8015003456163883`*^9}, {
   3.80150042807729*^9, 3.8015004684214773`*^9}, {3.801500550208126*^9, 
   3.8015005680612087`*^9}, {3.8015006719995937`*^9, 3.801500699909223*^9}, {
   3.8015007539132223`*^9, 3.801500807800034*^9}, {3.801500843269526*^9, 
   3.801500846056826*^9}, 3.801664244151732*^9, 3.801664286300025*^9, 
   3.801666771044989*^9, 3.801667066434883*^9, 3.8024505336003647`*^9},
 FontSize->12,
 CellLabel->"In[50]:=",ExpressionUUID->"ad49dce0-a7fe-4277-9a31-c57ad43d9c4e"],

Cell[BoxData[
 RowBox[{"a1", "=", 
  RowBox[{
   RowBox[{"Tan", "[", "\[Alpha]1", "]"}], "*", "g1", " ", 
   RowBox[{"(*", " ", 
    RowBox[{
     RowBox[{
      RowBox[{
      "then", " ", "m1", " ", "is", " ", "at", " ", "rest", " ", "in", " ", 
       "non"}], "-", 
      RowBox[{"inertial", " ", "frame"}]}], ";", " ", 
     RowBox[{
      RowBox[{
       RowBox[{"if", " ", "a1"}], "\[Equal]", 
       RowBox[{
       "0", " ", "then", " ", "the", " ", "answer", " ", "for", " ", "m1"}]}],
       "=", 
      RowBox[{
      "1", " ", "motion", " ", "must", " ", "be", " ", "as", " ", 
       "above"}]}]}], " ", "*)"}]}]}]], "Input",
 CellChangeTimes->{{3.8016641339809504`*^9, 3.801664178721177*^9}, {
   3.8016682988509703`*^9, 3.801668311594327*^9}, 3.801668394208872*^9, {
   3.801673909850916*^9, 3.801673917361527*^9}, {3.801674242150024*^9, 
   3.801674243741444*^9}, {3.8016742915471773`*^9, 3.8016742922672853`*^9}, {
   3.801674538820476*^9, 3.8016745412465878`*^9}, {3.8016746573244743`*^9, 
   3.801674667425908*^9}, {3.801674758299034*^9, 3.801674759522686*^9}, {
   3.801710500841217*^9, 3.801710550870118*^9}, {3.8017107211632423`*^9, 
   3.801710721608687*^9}, {3.8024505401960917`*^9, 3.802450557535974*^9}},
 CellLabel->"In[51]:=",ExpressionUUID->"4dcd7886-c79d-4389-bfcd-3d40fd8c3417"],

Cell[BoxData[
 RowBox[{"gM1", "=", 
  RowBox[{"Sqrt", "[", 
   RowBox[{
    RowBox[{"g1", "^", "2"}], "+", 
    RowBox[{"a1", "^", "2"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.801500135827752*^9, 3.801500156519579*^9}, {
  3.8016740685501633`*^9, 3.801674074886196*^9}, {3.801674139523796*^9, 
  3.801674150901403*^9}, {3.802450565259976*^9, 3.802450569607071*^9}},
 CellLabel->"In[52]:=",ExpressionUUID->"c4fb0e91-6165-42ff-bd0a-16a892ae8a00"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     RowBox[{"Tan", "[", "\[Beta]", "]"}], "=", 
     RowBox[{
      RowBox[{"Abs", "[", "a", "]"}], "/", "g"}]}], ";"}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"m", "=", "1"}], ";", "\n", 
    RowBox[{"\[Beta]", "=", 
     RowBox[{"ArcTan", "[", 
      RowBox[{
       RowBox[{"Abs", "[", "a", "]"}], "/", "g"}], "]"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"N", "[", "\[Beta]", "]"}], "\[IndentingNewLine]", 
      RowBox[{"RotationMatrix", "[", "\[Beta]", "]"}]}], "//", 
     RowBox[{"MatrixForm", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"RotationMatrix", "[", "\[Beta]", "]"}], ".", 
       RowBox[{"{", 
        RowBox[{"1", ",", "0"}], "}"}]}], "\[IndentingNewLine]", 
      RowBox[{"F", "=", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"m", "*", "g", "*", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"Cos", "[", 
             RowBox[{"\[Alpha]", "+", "\[Beta]"}], "]"}], "*", 
            RowBox[{"Sin", "[", 
             RowBox[{"\[Alpha]", "+", "\[Beta]"}], "]"}]}], ")"}]}], ",", 
         RowBox[{"m", "*", "g", "*", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{"Cos", "[", 
               RowBox[{"\[Alpha]", "+", "\[Beta]"}], "]"}], ")"}], "^", "2"}],
             "-", "1"}], ")"}]}]}], "}"}]}]}]}], ";", "\n", 
    RowBox[{"t1", "=", "0"}], ";", "\n", 
    RowBox[{"t2", "=", "3"}], ";", "\n", 
    RowBox[{"P0", "=", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"0", ",", 
        RowBox[{
         RowBox[{"-", "g"}], "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{
            RowBox[{"(", 
             RowBox[{"Cos", "[", "\[Alpha]", "]"}], ")"}], "^", "2"}], "-", 
           "1"}], ")"}]}]}], "}"}], "*", 
      RowBox[{
       RowBox[{"t2", "^", "2"}], "/", "2"}]}]}], ";", " ", 
    RowBox[{"(*", " ", 
     RowBox[{
     "adjusting", " ", "altitude", " ", "to", " ", "falling", " ", "to", " ", 
      "zero"}], " ", "*)"}], "\n", 
    RowBox[{"V0", "=", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0"}], "}"}]}], ";", "\n", 
    RowBox[{"steps", "=", "5000"}], ";", "\n", 
    RowBox[{"tf", "=", "0"}], ";", "\n", 
    RowBox[{"Motions", "[", 
     RowBox[{
     "m", ",", "unk", ",", "F", ",", "P0", ",", "V0", ",", "t1", ",", "t2", 
      ",", "steps", ",", "tf"}], "]"}]}], "\[IndentingNewLine]", 
   "*)"}]}]], "Input",
 CellChangeTimes->{{3.801664298581058*^9, 3.8016643065949287`*^9}, {
   3.8016643385334377`*^9, 3.801664408703248*^9}, {3.801664441359417*^9, 
   3.8016644427067537`*^9}, {3.8016644743550577`*^9, 3.801664606389697*^9}, 
   3.8016668109926033`*^9, {3.801666917429933*^9, 3.801666918092525*^9}, 
   3.801666955549739*^9, {3.801667073545068*^9, 3.801667094074245*^9}, {
   3.8016672282778788`*^9, 3.8016672500578337`*^9}, 3.801667433764526*^9, {
   3.80166752831491*^9, 3.801667534940559*^9}},
 CellLabel->"In[53]:=",ExpressionUUID->"8c2a8527-99b1-4999-8ec7-13a9dff68758"],

Cell[BoxData[
 RowBox[{"unk", "=", 
  RowBox[{"{", 
   RowBox[{"x", ",", "y"}], "}"}], " ", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"here", " ", "x", " ", "and", " ", "y"}], " ", "-", " ", 
    RowBox[{"usual", " ", "Cartesian"}]}], " ", "*)"}]}]], "Input",
 CellLabel->"In[54]:=",ExpressionUUID->"732e3690-6c47-4bc2-a3e3-399ec1db597b"],

Cell[BoxData[
 RowBox[{"\[Beta]1", "=", 
  RowBox[{"ArcTan", "[", 
   RowBox[{
    RowBox[{"-", "a1"}], "/", "g1"}], "]"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"minus", " ", "because", " ", "in", " ", "non"}], "-", 
    RowBox[{"inertial", " ", "frame", " ", "it", " ", "gives"}], " ", "-", 
    "a"}], " ", "*)"}]}]], "Input",
 CellChangeTimes->{{3.8016687028399153`*^9, 3.801668705012409*^9}, {
  3.801672626922172*^9, 3.801672673547915*^9}, {3.8024505784283943`*^9, 
  3.8024505833711042`*^9}},
 CellLabel->"In[55]:=",ExpressionUUID->"86949b1a-a94f-40fe-a008-106fcfaeae4d"],

Cell[BoxData[
 RowBox[{"t1", "=", "0"}]], "Input",
 CellChangeTimes->{{3.8016675941373663`*^9, 3.801667596768306*^9}},
 CellLabel->"In[56]:=",ExpressionUUID->"55d9fa88-0013-43f6-85eb-5f40a1c1baa0"],

Cell[BoxData[
 RowBox[{"t2", "=", "3"}]], "Input",
 CellChangeTimes->{{3.801667599595746*^9, 3.801667601560302*^9}},
 CellLabel->"In[57]:=",ExpressionUUID->"3d954a82-d59d-4213-ad29-836fb4c4a2b0"],

Cell[BoxData[
 RowBox[{"FnonIn1", "=", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"m1", "*", "gM1", "*", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"Cos", "[", 
        RowBox[{"\[Alpha]1", "+", "\[Beta]1"}], "]"}], "*", 
       RowBox[{"Sin", "[", 
        RowBox[{"\[Alpha]1", "+", "\[Beta]1"}], "]"}]}], ")"}]}], ",", 
    RowBox[{"m1", "*", "gM1", "*", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"Cos", "[", 
          RowBox[{"\[Alpha]1", "+", "\[Beta]1"}], "]"}], ")"}], "^", "2"}], 
       "-", "1"}], ")"}]}]}], "}"}]}]], "Input",
 CellChangeTimes->{{3.8016679580057898`*^9, 3.801667962446694*^9}, {
  3.801673771557877*^9, 3.8016737849756937`*^9}, {3.802450590694397*^9, 
  3.802450627360358*^9}},
 CellLabel->"In[58]:=",ExpressionUUID->"c2d1a6fc-b0e3-49c3-b5af-9d3ed52b0f86"],

Cell[BoxData[
 RowBox[{"Na1", "=", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"m1", "*", "gM1", "*", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"Cos", "[", 
        RowBox[{"\[Alpha]1", "+", "\[Beta]1"}], "]"}], "*", 
       RowBox[{"Sin", "[", 
        RowBox[{"\[Alpha]1", "+", "\[Beta]1"}], "]"}]}], ")"}]}], ",", 
    RowBox[{"m1", "*", "gM1", "*", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"Cos", "[", 
        RowBox[{"\[Alpha]1", "+", "\[Beta]1"}], "]"}], ")"}], "^", "2"}]}]}], 
   "}"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{
   "\:0440\:0435\:0430\:043a\:0446\:0438\:044f", " ", 
    "\:043e\:043f\:043e\:0440\:044b", " ", "\:0432", " ", 
    "\:0441\:0438\:0441\:0442\:0435\:043c\:0435", " ", "a"}], " ", 
   "*)"}]}]], "Input",
 CellChangeTimes->{{3.8016679580057898`*^9, 3.801667962446694*^9}, {
  3.801668045829935*^9, 3.801668118098494*^9}, {3.802450630523918*^9, 
  3.802450649178774*^9}},
 CellLabel->"In[59]:=",ExpressionUUID->"513f85a1-1933-4b28-a727-4c13f8e32bc1"],

Cell[BoxData[
 RowBox[{
  RowBox[{"RotationMatrix", "[", "\[Beta]1", "]"}], "//", 
  "MatrixForm"}]], "Input",
 CellChangeTimes->{3.8024506540032043`*^9},
 CellLabel->"In[60]:=",ExpressionUUID->"aabe84ee-187b-4050-aaae-51d1d8188fad"],

Cell[BoxData[
 RowBox[{"N01", "=", 
  RowBox[{
   RowBox[{"RotationMatrix", "[", "\[Beta]1", "]"}], ".", "Na1", " ", 
   RowBox[{"(*", " ", 
    RowBox[{
    "\:0440\:0435\:0430\:043a\:0446\:0438\:044f", " ", 
     "\:043e\:043f\:043e\:0440\:044b", " ", "\:0432", " ", 
     "\:0438\:0441\:0445\:043e\:0434\:043d\:044b\:0445", " ", 
     "\:043a\:043e\:043e\:0440\:0434\:0438\:043d\:0430\:0442\:0430\:0445"}], 
    " ", "*)"}]}]}]], "Input",
 CellChangeTimes->{{3.801668210758245*^9, 3.80166822643504*^9}, {
  3.801668259694399*^9, 3.8016682855321293`*^9}, {3.802450662587696*^9, 
  3.802450668000408*^9}},
 CellLabel->"In[61]:=",ExpressionUUID->"1d3a1415-ca3a-4501-8722-26383d7543ed"],

Cell[BoxData[
 RowBox[{"F1", "=", " ", 
  RowBox[{"N01", "+", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{"-", "m1"}], "*", "a1"}], ",", 
     RowBox[{
      RowBox[{"-", "m1"}], "*", "g1"}]}], "}"}]}]}]], "Input",
 CellChangeTimes->{{3.8016679580057898`*^9, 3.801667962446694*^9}, {
  3.801673677174209*^9, 3.801673754568101*^9}, {3.802450680285242*^9, 
  3.80245068540213*^9}},
 CellLabel->"In[62]:=",ExpressionUUID->"c1fbf734-e9ce-4f40-8932-6aadfcdcc4de"],

Cell[BoxData[
 RowBox[{"P01", "=", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"0", ",", 
     RowBox[{
      RowBox[{"-", "g1"}], "*", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{"Cos", "[", "\[Alpha]1", "]"}], ")"}], "^", "2"}], "-", 
        "1"}], ")"}]}]}], "}"}], "*", 
   RowBox[{
    RowBox[{"t2", "^", "2"}], "/", "2"}]}]}]], "Input",
 CellChangeTimes->{{3.802450693638896*^9, 3.802450698262323*^9}},
 CellLabel->"In[63]:=",ExpressionUUID->"f06ae692-8295-4c7e-905a-e4d5cea4a4aa"],

Cell[BoxData[
 RowBox[{"V01", "=", 
  RowBox[{"{", 
   RowBox[{"0", ",", "0"}], "}"}]}]], "Input",
 CellChangeTimes->{3.802450708545555*^9},
 CellLabel->"In[64]:=",ExpressionUUID->"14472c62-f4dc-4610-b5bd-a934de995d72"],

Cell[BoxData[
 RowBox[{"steps", "=", "5000"}]], "Input",
 CellChangeTimes->{{3.801667745754609*^9, 3.801667749197032*^9}},
 CellLabel->"In[65]:=",ExpressionUUID->"9d7133bb-c1d5-4c68-b218-809bad881501"],

Cell[BoxData[
 RowBox[{"tf", "=", "0"}]], "Input",
 CellChangeTimes->{{3.801667762253001*^9, 3.801667763747047*^9}},
 CellLabel->"In[66]:=",ExpressionUUID->"84e4be1e-158f-4191-acf5-e481203f979e"],

Cell[BoxData[
 RowBox[{"Motions", "[", 
  RowBox[{
  "m1", ",", "unk", ",", "F1", ",", "P01", ",", "V01", ",", "t1", ",", "t2", 
   ",", "steps", ",", "tf"}], "]"}]], "Input",
 CellChangeTimes->{{3.802450715045269*^9, 3.802450718810597*^9}},
 CellLabel->"In[67]:=",ExpressionUUID->"9f9d6631-b1cb-4326-bfa6-c7d732cf2f38"]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
 "E",
 StyleBox["xample",
  FontFamily->"Times New Roman"],
 " seb4: ",
 StyleBox["Material point with mass ",
  FontWeight->"Bold"],
 StyleBox["m",
  FontFamily->"Utopia",
  FontSlant->"Italic"],
 StyleBox[" moving along an inclined (angle \[Alpha]) plain (wedge) on \
straight line of steepest gradient and acted upon by gravity and reaction of \
surface. The wedge of mass ",
  FontWeight->"Bold"],
 StyleBox["M",
  FontFamily->"Utopia",
  FontSlant->"Italic"],
 StyleBox[" is moving without friction on table. Find acceleration ",
  FontWeight->"Bold"],
 StyleBox["a",
  FontFamily->"Utopia",
  FontSlant->"Italic"],
 StyleBox[" of the wedge (",
  FontWeight->"Bold"],
 StyleBox["\:041c\:0435\:0449\:0435\:0440\:0441\:043a\:0438\:0439",
  FontFamily->"Ubuntu Mono",
  FontWeight->"Bold"],
 StyleBox[" 48.28\.13) . ",
  FontWeight->"Bold"]
}], "Subsection",
 Evaluatable->False,
 CellChangeTimes->{{3.7227623633962064`*^9, 3.7227623643477526`*^9}, {
   3.7227624078391485`*^9, 3.722762412550025*^9}, 3.7239814420362844`*^9, {
   3.8014125894411163`*^9, 3.80141259052518*^9}, {3.801412731179865*^9, 
   3.801412742256505*^9}, 3.80149934928935*^9, {3.8014994372951393`*^9, 
   3.801499559283052*^9}, {3.801499610578313*^9, 3.801499629038383*^9}, {
   3.801499940257044*^9, 3.8014999526139793`*^9}, {3.801666617251504*^9, 
   3.801666618525972*^9}, {3.8021495159595203`*^9, 3.8021495180469513`*^9}, {
   3.8021495575449057`*^9, 3.802149652872552*^9}, {3.802149840446677*^9, 
   3.8021498631530943`*^9}},
 AspectRatioFixed->True,
 FontFamily->"Times New Roman",
 FontSize->16,ExpressionUUID->"d92c15d0-00fb-4ff4-aabb-3632c6ddb818"],

Cell["\<\
Consider ref. frame where wedge is at rest. Then it is non-inertial and we \
have equiv. gravity \
\>", "Text",
 CellChangeTimes->{{3.801500034245212*^9, 
  3.801500093033792*^9}},ExpressionUUID->"2868ac19-adf4-4edc-92df-\
e570bddf8523"],

Cell[TextData[{
 "Find acceleration ",
 StyleBox["a1",
  FontSlant->"Italic"],
 " for",
 StyleBox[" M",
  FontSlant->"Italic"],
 " when ",
 StyleBox["m",
  FontSlant->"Italic"],
 " is sliding down."
}], "Text",
 CellChangeTimes->{{3.801674944265684*^9, 3.801674994361802*^9}, {
  3.802104331604205*^9, 3.802104369156797*^9}, {3.8021501087785482`*^9, 
  3.802150114254904*^9}, {3.802449505065085*^9, 3.802449529242771*^9}, {
  3.8024499322606564`*^9, 
  3.802449937007093*^9}},ExpressionUUID->"e3b11540-24d8-450f-9c78-\
55724d87f143"],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{"\[Alpha]", "=", 
   RowBox[{"\[Pi]", "/", "6"}]}], " ", "*)"}]], "Input",
 CellChangeTimes->{{3.488516603965781*^9, 3.488516609450156*^9}, 
   3.533638409953125*^9, {3.533639196609375*^9, 3.533639200109375*^9}, 
   3.533639232609375*^9, 3.7227627607774453`*^9, {3.7239814547375064`*^9, 
   3.723981456019337*^9}, {3.8014127751361628`*^9, 3.8014128097557993`*^9}, {
   3.801412840033223*^9, 3.801412840844523*^9}, {3.8014129630674887`*^9, 
   3.801412974090295*^9}, 3.801413298214366*^9, {3.801413365263443*^9, 
   3.8014133711072283`*^9}, {3.801413500442295*^9, 3.801413521750505*^9}, {
   3.801413766388073*^9, 3.8014138548879213`*^9}, {3.80141397735352*^9, 
   3.801413996378757*^9}, {3.8014983672015257`*^9, 3.801498394810019*^9}, {
   3.801498484419816*^9, 3.801498516950513*^9}, {3.801498550119773*^9, 
   3.8014985530588493`*^9}, {3.801498639396269*^9, 3.8014986529798727`*^9}, {
   3.801498684739962*^9, 3.8014987056069527`*^9}, {3.8014987616399307`*^9, 
   3.801498840921955*^9}, 3.801498880172531*^9, {3.801498920206132*^9, 
   3.801498988748966*^9}, {3.801499084421528*^9, 3.801499084634947*^9}, {
   3.801499576164363*^9, 3.8014995763877296`*^9}, {3.801500116322962*^9, 
   3.8015001235638103`*^9}, {3.801500297090946*^9, 3.8015003456163883`*^9}, {
   3.80150042807729*^9, 3.8015004684214773`*^9}, {3.801500550208126*^9, 
   3.8015005680612087`*^9}, {3.8015006719995937`*^9, 3.801500699909223*^9}, {
   3.8015007539132223`*^9, 3.801500807800034*^9}, {3.801500843269526*^9, 
   3.801500846056826*^9}, 3.801664244151732*^9, 3.801664286300025*^9, 
   3.801666771044989*^9, 3.801667066434883*^9, {3.802449824383357*^9, 
   3.802449828506946*^9}},
 FontSize->12,
 CellLabel->"In[68]:=",ExpressionUUID->"6731fb10-3492-4361-bae2-795baac49b00"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"a", "=", 
    RowBox[{
     RowBox[{"Tan", "[", "\[Alpha]", "]"}], "*", "g"}]}], " ", "*)"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{
     RowBox[{
     "then", " ", "m", " ", "is", " ", "at", " ", "rest", " ", "in", " ", 
      "non"}], "-", 
     RowBox[{"inertial", " ", "frame"}]}], ";", " ", 
    RowBox[{
     RowBox[{
      RowBox[{"if", " ", "a"}], "\[Equal]", 
      RowBox[{
      "0", " ", "then", " ", "the", " ", "answer", " ", "for", " ", "m"}]}], 
     "=", 
     RowBox[{
     "1", " ", "motion", " ", "must", " ", "be", " ", "as", " ", 
      "above"}]}]}], " ", "*)"}]}]], "Input",
 CellChangeTimes->{{3.8016641339809504`*^9, 3.801664178721177*^9}, {
   3.8016682988509703`*^9, 3.801668311594327*^9}, 3.801668394208872*^9, {
   3.801673909850916*^9, 3.801673917361527*^9}, {3.801674242150024*^9, 
   3.801674243741444*^9}, {3.8016742915471773`*^9, 3.8016742922672853`*^9}, {
   3.801674538820476*^9, 3.8016745412465878`*^9}, {3.8016746573244743`*^9, 
   3.801674667425908*^9}, {3.801674758299034*^9, 3.801674759522686*^9}, {
   3.801710500841217*^9, 3.801710550870118*^9}, {3.8017107211632423`*^9, 
   3.801710721608687*^9}, {3.80244990058278*^9, 3.8024499153336906`*^9}},
 CellLabel->"In[69]:=",ExpressionUUID->"8ec0362c-72d0-430d-bb9f-26e15a26d490"],

Cell[BoxData[
 RowBox[{"gM", "=", 
  RowBox[{"Sqrt", "[", 
   RowBox[{
    RowBox[{"g", "^", "2"}], "+", 
    RowBox[{"a", "^", "2"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.801500135827752*^9, 3.801500156519579*^9}, {
   3.8016740685501633`*^9, 3.801674074886196*^9}, {3.801674139523796*^9, 
   3.801674150901403*^9}, 3.802449943543969*^9, 3.8024500743316727`*^9, {
   3.802451224059806*^9, 3.8024512266470613`*^9}},
 CellLabel->"In[70]:=",ExpressionUUID->"1fc4292c-3ccd-4980-a1db-09e0e4ea8d77"],

Cell[BoxData[
 RowBox[{"unk", "=", 
  RowBox[{"{", 
   RowBox[{"x", ",", "y"}], "}"}], " ", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"here", " ", "x", " ", "and", " ", "y"}], " ", "-", " ", 
    RowBox[{"usual", " ", "Cartesian"}]}], " ", "*)"}]}]], "Input",
 CellLabel->"In[71]:=",ExpressionUUID->"63a0e785-9b67-497d-b91a-0bd02737f2ac"],

Cell[BoxData[
 RowBox[{"\[Beta]", "=", 
  RowBox[{"ArcTan", "[", 
   RowBox[{
    RowBox[{"-", "a"}], "/", "g"}], "]"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"minus", " ", "because", " ", "in", " ", "non"}], "-", 
    RowBox[{"inertial", " ", "frame", " ", "it", " ", "gives"}], " ", "-", 
    "a"}], " ", "*)"}]}]], "Input",
 CellChangeTimes->{{3.8016687028399153`*^9, 3.801668705012409*^9}, {
   3.801672626922172*^9, 3.801672673547915*^9}, 3.8024499491596403`*^9, 
   3.802450080993197*^9, {3.8024512326760073`*^9, 3.802451234373794*^9}},
 CellLabel->"In[72]:=",ExpressionUUID->"3ac44665-b330-411c-8301-435d36908828"],

Cell[BoxData[
 RowBox[{"t1", "=", "0"}]], "Input",
 CellChangeTimes->{{3.8016675941373663`*^9, 3.801667596768306*^9}},
 CellLabel->"In[73]:=",ExpressionUUID->"f7cf8411-0145-4367-9805-85220d4a15e9"],

Cell[BoxData[
 RowBox[{"t2", "=", "3"}]], "Input",
 CellChangeTimes->{{3.801667599595746*^9, 3.801667601560302*^9}},
 CellLabel->"In[74]:=",ExpressionUUID->"bce37c15-84c9-4591-93a8-54fd7445d104"],

Cell[TextData[{
 "Force acting on ",
 StyleBox["m",
  FontSlant->"Italic"],
 " in non-inertial frame:"
}], "Text",
 CellChangeTimes->{{3.802449666597499*^9, 
  3.8024496966590767`*^9}},ExpressionUUID->"2c9309f8-1597-4f85-8bef-\
5d9f88aacc04"],

Cell[BoxData[
 RowBox[{"FnonIn", "=", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"m", "*", "gM", "*", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"Cos", "[", 
        RowBox[{"\[Alpha]", "+", "\[Beta]"}], "]"}], "*", 
       RowBox[{"Sin", "[", 
        RowBox[{"\[Alpha]", "+", "\[Beta]"}], "]"}]}], ")"}]}], ",", 
    RowBox[{"m", "*", "gM", "*", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"Cos", "[", 
          RowBox[{"\[Alpha]", "+", "\[Beta]"}], "]"}], ")"}], "^", "2"}], "-",
        "1"}], ")"}]}]}], "}"}]}]], "Input",
 CellChangeTimes->{{3.8016679580057898`*^9, 3.801667962446694*^9}, {
  3.801673771557877*^9, 3.8016737849756937`*^9}},
 CellLabel->"In[75]:=",ExpressionUUID->"0a116422-d0d9-41a9-bdc2-f938ef32d89a"],

Cell[BoxData[
 RowBox[{"Na", "=", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"m", "*", "gM", "*", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"Cos", "[", 
        RowBox[{"\[Alpha]", "+", "\[Beta]"}], "]"}], "*", 
       RowBox[{"Sin", "[", 
        RowBox[{"\[Alpha]", "+", "\[Beta]"}], "]"}]}], ")"}]}], ",", 
    RowBox[{"m", "*", "gM", "*", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"Cos", "[", 
        RowBox[{"\[Alpha]", "+", "\[Beta]"}], "]"}], ")"}], "^", "2"}]}]}], 
   "}"}], " ", 
  RowBox[{"(*", " ", 
   StyleBox[
    RowBox[{
    "\:0440\:0435\:0430\:043a\:0446\:0438\:044f", " ", 
     "\:043e\:043f\:043e\:0440\:044b", " ", "\:0432", " ", 
     "\:0441\:0438\:0441\:0442\:0435\:043c\:0435", " ", 
     "\:043f\:043e\:043a\:043e\:044f", " ", "\:043c\:0430\:0441\:0441\:044b", 
     " ", "M"}],
    FontFamily->"Ubuntu Mono",
    FontWeight->"Plain"], " ", "*)"}]}]], "Input",
 CellChangeTimes->{{3.8016679580057898`*^9, 3.801667962446694*^9}, {
  3.801668045829935*^9, 3.801668118098494*^9}, {3.8023332684569674`*^9, 
  3.802333270900882*^9}, {3.802333302648278*^9, 3.80233330910367*^9}},
 CellLabel->"In[76]:=",ExpressionUUID->"1617cc33-a1a6-4c80-bf98-d6a813d7eb14"],

Cell[BoxData[
 RowBox[{
  RowBox[{"RotationMatrix", "[", "\[Beta]", "]"}], "//", 
  "MatrixForm"}]], "Input",
 CellLabel->"In[77]:=",ExpressionUUID->"59f083c4-7601-459c-b5b3-93e7072b5f2f"],

Cell[BoxData[
 RowBox[{"N0", "=", 
  RowBox[{
   RowBox[{"RotationMatrix", "[", "\[Beta]", "]"}], ".", "Na", " ", 
   RowBox[{"(*", " ", 
    RowBox[{
    "\:0440\:0435\:0430\:043a\:0446\:0438\:044f", " ", 
     "\:043e\:043f\:043e\:0440\:044b", " ", "\:0432", " ", 
     "\:0438\:0441\:0445\:043e\:0434\:043d\:044b\:0445", " ", 
     "\:043a\:043e\:043e\:0440\:0434\:0438\:043d\:0430\:0442\:0430\:0445"}], 
    " ", "*)"}]}]}]], "Input",
 CellChangeTimes->{{3.801668210758245*^9, 3.80166822643504*^9}, {
  3.801668259694399*^9, 3.8016682855321293`*^9}},
 CellLabel->"In[78]:=",ExpressionUUID->"1c49f1f1-c13e-4999-9be6-e7336c8513c2"],

Cell[BoxData[
 RowBox[{"F", "=", " ", 
  RowBox[{"N0", "+", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{"-", "m"}], "*", "a"}], ",", 
     RowBox[{
      RowBox[{"-", "m"}], "*", "g"}]}], "}"}]}]}]], "Input",
 CellChangeTimes->{{3.8016679580057898`*^9, 3.801667962446694*^9}, {
  3.801673677174209*^9, 3.801673754568101*^9}, {3.802450116910671*^9, 
  3.8024501192693243`*^9}, {3.802451253461335*^9, 3.802451255492812*^9}},
 CellLabel->"In[79]:=",ExpressionUUID->"b6e9f689-1e7a-4b7a-9e57-79620e8bf89e"],

Cell[BoxData[
 RowBox[{"FullSimplify", "[", "F", "]"}]], "Input",
 CellChangeTimes->{{3.802451452672119*^9, 3.8024514628584948`*^9}, {
  3.802451515445614*^9, 3.802451516593781*^9}},
 CellLabel->"In[80]:=",ExpressionUUID->"ed5fd72c-853d-424a-91b8-99f5d2a46e1c"],

Cell["\<\
Hence, in original reference frame we have acceleration of mass m:\
\>", "Text",
 CellChangeTimes->{{3.8024992194243097`*^9, 
  3.802499257165635*^9}},ExpressionUUID->"fcded526-a732-47b0-a93d-\
3bc366341383"],

Cell[BoxData[
 RowBox[{"a_m", "=", " ", 
  FractionBox[
   RowBox[{
    SqrtBox[
     RowBox[{"1", "+", 
      FractionBox[
       SuperscriptBox["a", "2"], 
       SuperscriptBox["g", "2"]]}]], " ", "g", "  ", 
    RowBox[{"Sin", "[", "\[Alpha]", "]"}], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"g", " ", 
       RowBox[{"Cos", "[", "\[Alpha]", "]"}]}], "+", 
      RowBox[{"a", " ", 
       RowBox[{"Sin", "[", "\[Alpha]", "]"}]}]}], ")"}]}], 
   SqrtBox[
    RowBox[{
     SuperscriptBox["a", "2"], "+", 
     SuperscriptBox["g", "2"]}]]]}]], "Input",
 CellChangeTimes->{{3.802499269024618*^9, 3.802499307863164*^9}},
 CellLabel->"In[81]:=",ExpressionUUID->"c629bbb1-11a0-488c-8833-6bb096cc68aa"],

Cell["But due to the conservation of momentum we must have", "Text",
 CellChangeTimes->{{3.802499387784659*^9, 3.80249940703376*^9}, 
   3.8024994922356253`*^9},ExpressionUUID->"664e13cf-25b1-448c-bf9f-\
83c8c916553b"],

Cell[BoxData[
 RowBox[{"Solve", "[", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"m", "*", 
      FractionBox[
       RowBox[{
        SqrtBox[
         RowBox[{"1", "+", 
          FractionBox[
           SuperscriptBox["a", "2"], 
           SuperscriptBox["g", "2"]]}]], " ", "g", "  ", 
        RowBox[{"Sin", "[", "\[Alpha]", "]"}], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"g", " ", 
           RowBox[{"Cos", "[", "\[Alpha]", "]"}]}], "+", 
          RowBox[{"a", " ", 
           RowBox[{"Sin", "[", "\[Alpha]", "]"}]}]}], ")"}]}], 
       SqrtBox[
        RowBox[{
         SuperscriptBox["a", "2"], "+", 
         SuperscriptBox["g", "2"]}]]]}], "+", 
     RowBox[{"M", "*", "a"}]}], "==", "0"}], ",", "a"}], "]"}]], "Input",
 CellChangeTimes->{{3.8024994466431913`*^9, 3.80249945924886*^9}, {
  3.802499508074991*^9, 3.802499524342792*^9}, {3.802499559777329*^9, 
  3.802499565369158*^9}, {3.802499615449627*^9, 3.8024996218823338`*^9}},
 CellLabel->"In[82]:=",ExpressionUUID->"6001bdcc-5a9a-4b01-91a9-db19fb9171d8"],

Cell["\<\
From \[OpenCurlyDoubleQuote]Introduction to Classical Mechanics With Problems \
and Solutions\[CloseCurlyDoubleQuote]
David Morin, Cambridge University Press - 2007
 eq. (3.83)

Ax = mg sin \[Theta] cos \[Theta]/(M + m sin^2 \[Theta])
 
 \:0410. \:0412. \:041f\:0430\:043d\:0448\:0438\:043d\:0430, \:0412. \:041c. \
\:0427\:0443\:0440\:043a\:0438\:043d
 \:0422\:0415\:041e\:0420\:0415\:0422\:0418\:0427\:0415\:0421\:041a\:0410\
\:042f \:041c\:0415\:0425\:0410\:041d\:0418\:041a\:0410 \:0412 \:0420\:0415\
\:0428\:0415\:041d\:0418\:042f\:0425 \:0417\:0410\:0414\:0410\:0427 \:0418\
\:0417 \:0421\:0411\:041e\:0420\:041d\:0418\:041a\:0410 \:0418. \:0412. \
\:041c\:0415\:0429\:0415\:0420\:0421\:041a\:041e\:0413\:041e
 \:0421\:0442\:0440. 141. \:0417\:0430\:0434\:0430\:0447\:0430 3.22(48.28).
 mg sin 2\[Alpha] /[ 2 (m_1 +m sin^2 \[Alpha])]
\
\>", "Text",
 CellChangeTimes->{{3.802499888233379*^9, 3.802499904487105*^9}, {
  3.802499936535244*^9, 3.8024999834426403`*^9}, {3.802500027891179*^9, 
  3.8025000880808487`*^9}, {3.8025001747482367`*^9, 3.802500203262632*^9}, {
  3.80250023808528*^9, 3.802500285791088*^9}, {3.802500652625589*^9, 
  3.802500695174492*^9}, {3.802503987758606*^9, 3.802504001307959*^9}, {
  3.802504041615613*^9, 3.802504059460782*^9}, {3.802505235177577*^9, 
  3.802505268133882*^9}},ExpressionUUID->"efc5cf6b-ff48-4f83-833c-\
24a863cac639"],

Cell[BoxData[
 RowBox[{"FullSimplify", "[", 
  FractionBox[
   RowBox[{"g", " ", "m", " ", 
    RowBox[{"Sin", "[", 
     RowBox[{"2", " ", "\[Alpha]"}], "]"}]}], 
   RowBox[{
    RowBox[{"-", "m"}], "-", 
    RowBox[{"2", " ", "M"}], "+", 
    RowBox[{"m", " ", 
     RowBox[{"Cos", "[", 
      RowBox[{"2", " ", "\[Alpha]"}], "]"}]}]}]], "]"}]], "Input",
 CellChangeTimes->{{3.802500702351305*^9, 3.802500728180819*^9}},
 CellLabel->"In[83]:=",ExpressionUUID->"a270da95-89f6-4280-b600-e6a07cf58d65"],

Cell[BoxData[
 RowBox[{"TrigExpand", "[", 
  FractionBox[
   RowBox[{"g", " ", "m", " ", 
    RowBox[{"Sin", "[", 
     RowBox[{"2", " ", "\[Alpha]"}], "]"}]}], 
   RowBox[{
    RowBox[{"-", "m"}], "-", 
    RowBox[{"2", " ", "M"}], "+", 
    RowBox[{"m", " ", 
     RowBox[{"Cos", "[", 
      RowBox[{"2", " ", "\[Alpha]"}], "]"}]}]}]], "]"}]], "Input",
 NumberMarks->False,
 CellLabel->"In[84]:=",ExpressionUUID->"07ada612-cb0d-4a82-9cbb-d9fb2b98a5af"],

Cell[TextData[{
 Cell[BoxData[
  FractionBox[
   RowBox[{"2", " ", "g", " ", "m", " ", 
    RowBox[{"Cos", "[", "\[Alpha]", "]"}], " ", 
    RowBox[{"Sin", "[", "\[Alpha]", "]"}]}], 
   RowBox[{
    RowBox[{"-", "m"}], "-", 
    RowBox[{"2", " ", "M"}], "+", 
    RowBox[{"m", " ", 
     SuperscriptBox[
      RowBox[{"Cos", "[", "\[Alpha]", "]"}], "2"]}], "-", 
    RowBox[{"m", " ", 
     SuperscriptBox[
      RowBox[{"Sin", "[", "\[Alpha]", "]"}], "2"]}]}]]], "Output",
  CellChangeTimes->{3.802503948422488*^9},ExpressionUUID->
  "934335ba-c0bf-44c7-b0d0-dedb6b4a0384"],
 "=",
 Cell[BoxData[
  FractionBox[
   RowBox[{"2", " ", "g", " ", "m", " ", 
    RowBox[{"Cos", "[", "\[Alpha]", "]"}], " ", 
    RowBox[{"Sin", "[", "\[Alpha]", "]"}]}], 
   RowBox[{
    RowBox[{
     RowBox[{"-", "2"}], " ", "M"}], " ", "-", 
    RowBox[{"m", 
     RowBox[{"(", 
      RowBox[{"1", "-", 
       SuperscriptBox[
        RowBox[{"Cos", "[", "\[Alpha]", "]"}], "2"]}], ")"}]}], "-", 
    RowBox[{"m", " ", 
     SuperscriptBox[
      RowBox[{"Sin", "[", "\[Alpha]", "]"}], "2"]}]}]]], "Output",
  GeneratedCell->False,
  CellAutoOverwrite->False,
  CellChangeTimes->{3.802503948422488*^9},ExpressionUUID->
  "7b42976f-1ea5-43c1-8a44-d49d520d1045"],
 "=",
 Cell[BoxData[
  FractionBox[
   RowBox[{"2", " ", "g", " ", "m", " ", 
    RowBox[{"Cos", "[", "\[Alpha]", "]"}], " ", 
    RowBox[{"Sin", "[", "\[Alpha]", "]"}]}], 
   RowBox[{
    RowBox[{
     RowBox[{"-", "2"}], " ", "M"}], " ", "-", 
    RowBox[{"2", "m", " ", 
     SuperscriptBox[
      RowBox[{"Sin", "[", "\[Alpha]", "]"}], "2"]}]}]]], "Output",
  GeneratedCell->False,
  CellAutoOverwrite->False,
  CellChangeTimes->{3.802503948422488*^9},ExpressionUUID->
  "36bf4690-ce44-4919-a1e7-1de334d95fec"],
 "=",
 Cell[BoxData[
  FractionBox[
   RowBox[{
    RowBox[{"-", " ", "g"}], " ", "m", " ", 
    RowBox[{"Sin", "[", 
     RowBox[{"2", "\[Alpha]"}], "]"}]}], 
   RowBox[{"2", 
    RowBox[{"(", " ", 
     RowBox[{"M", " ", "+", 
      RowBox[{"m", " ", 
       SuperscriptBox[
        RowBox[{"Sin", "[", "\[Alpha]", "]"}], "2"]}]}], ")"}]}]]], "Output",
  GeneratedCell->False,
  CellAutoOverwrite->False,
  CellChangeTimes->{3.802503948422488*^9},ExpressionUUID->
  "0e5bb8df-7ac9-4288-b99d-94f3a5fdf55a"],
 "\n\:0442.\:0435. \:043e\:0442\:0432\:0435\:0442\:044b \:0441\:043e\:0432\
\:043f\:0430\:0434\:0430\:044e\:0442!"
}], "Text",
 CellChangeTimes->{{3.802505430596395*^9, 
  3.802505621825457*^9}},ExpressionUUID->"6935b142-1c1c-464b-9871-\
f12ac94de75d"],

Cell[TextData[{
 "P0 = initial position of mass ",
 StyleBox["m",
  FontSlant->"Italic"],
 ":"
}], "Text",
 CellChangeTimes->{{3.8023334878497267`*^9, 
  3.8023335195199842`*^9}},ExpressionUUID->"b66c14d1-2065-4d85-acac-\
4a5dac7a8f6b"],

Cell[BoxData[
 RowBox[{"P0", "=", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"0", ",", 
     RowBox[{
      RowBox[{"-", "g"}], "*", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{"Cos", "[", "\[Alpha]", "]"}], ")"}], "^", "2"}], "-", 
        "1"}], ")"}]}]}], "}"}], "*", 
   RowBox[{
    RowBox[{"t2", "^", "2"}], "/", "2"}]}]}]], "Input",
 CellChangeTimes->{3.8024501297918863`*^9, 3.802451261838747*^9},
 CellLabel->"In[85]:=",ExpressionUUID->"4b636c6e-465d-4fe8-8d45-ed4f40e8522c"],

Cell[BoxData[
 RowBox[{"V0", "=", 
  RowBox[{"{", 
   RowBox[{"0", ",", "0"}], "}"}]}]], "Input",
 CellLabel->"In[86]:=",ExpressionUUID->"1986ff60-91c1-4b18-b28f-508c409cc8c6"],

Cell[BoxData[
 RowBox[{"steps", "=", "5000"}]], "Input",
 CellChangeTimes->{{3.801667745754609*^9, 3.801667749197032*^9}},
 CellLabel->"In[87]:=",ExpressionUUID->"1a6ae233-3ef8-404f-acbf-78216b4b1db1"],

Cell[BoxData[
 RowBox[{"tf", "=", "0"}]], "Input",
 CellChangeTimes->{{3.801667762253001*^9, 3.801667763747047*^9}},
 CellLabel->"In[88]:=",ExpressionUUID->"92f91791-ba56-469f-b7d6-5dd728ceb56e"],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{"Motions", "[", 
   RowBox[{
   "m", ",", "unk", ",", "F", ",", "P0", ",", "V0", ",", "t1", ",", "t2", ",",
     "steps", ",", "tf"}], "]"}], " ", "*)"}]], "Input",
 CellChangeTimes->{{3.802499061064116*^9, 3.802499065289221*^9}},
 CellLabel->"In[89]:=",ExpressionUUID->"bea3a17f-9df5-4f0b-9b42-6637438532df"]
}, Open  ]]
}, Open  ]]
}, Open  ]]
},
WindowSize->{1600, 829},
WindowMargins->{{0, Automatic}, {0, Automatic}},
Magnification->1.5,
FrontEndVersion->"11.3 for Linux x86 (64-bit) (March 6, 2018)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 544, 17, 231, "Subtitle",ExpressionUUID->"df9343b8-c2ab-43c3-8555-2ada59357c0d",
 Evaluatable->False],
Cell[1127, 41, 1804, 29, 190, "Text",ExpressionUUID->"9e441257-d59f-4c27-9fa8-64ea2349ba29"],
Cell[CellGroupData[{
Cell[2956, 74, 411, 12, 80, "Section",ExpressionUUID->"2e9c14ac-be8d-4a96-9d34-cc3f155d6741"],
Cell[3370, 88, 1052, 20, 47, "Text",ExpressionUUID->"7c93fce4-aeec-4fcd-8316-27364d5517c4"]
}, Open  ]],
Cell[CellGroupData[{
Cell[4459, 113, 422, 12, 80, "Section",ExpressionUUID->"092d8b89-2f67-40c9-a994-67dc31fd59a8"],
Cell[4884, 127, 728, 10, 47, "Text",ExpressionUUID->"04780f57-f191-4331-8c7c-9798e27441d0"],
Cell[5615, 139, 586, 15, 36, "DisplayFormula",ExpressionUUID->"844f2f14-e18e-4873-887a-4bcb4ce7f67b"],
Cell[6204, 156, 650, 10, 47, "Text",ExpressionUUID->"88a7d7c6-3f2f-4a34-8a1b-10caf5a44030"],
Cell[CellGroupData[{
Cell[6879, 170, 296, 4, 42, "Item",ExpressionUUID->"96ebc100-b7a2-4c58-9483-a751106aca3e"],
Cell[7178, 176, 440, 8, 42, "Item",ExpressionUUID->"ebd1d692-89e7-4ff8-8091-e01e2f1ce9bd"],
Cell[7621, 186, 299, 4, 42, "Item",ExpressionUUID->"0428a09f-be09-4e1c-8ffd-0c577684066b"],
Cell[7923, 192, 304, 4, 42, "Item",ExpressionUUID->"800b51e4-5a10-4045-ab95-b437f2cd7467"],
Cell[8230, 198, 304, 4, 42, "Item",ExpressionUUID->"48f54f2a-2bbe-4ae4-a715-924236549001"],
Cell[8537, 204, 311, 4, 42, "Item",ExpressionUUID->"36b72aa8-e58d-4a61-a4d9-e4760dfef2a2"],
Cell[8851, 210, 267, 3, 42, "Item",ExpressionUUID->"411034f7-dc86-44ef-8bf4-6c79bdd98508"],
Cell[9121, 215, 645, 15, 42, "Item",ExpressionUUID->"2823ec98-62bb-4cc1-88a4-c559fbecac7a"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[9815, 236, 428, 12, 80, "Section",ExpressionUUID->"be9ad61f-e973-403d-84b4-e561c54b767e"],
Cell[10246, 250, 13072, 324, 1453, "Input",ExpressionUUID->"626c1585-ac06-43dc-bf9d-f5a751a84769"]
}, Open  ]],
Cell[CellGroupData[{
Cell[23355, 579, 156, 2, 79, "Section",ExpressionUUID->"9d1eb348-1a14-4b6a-853a-7747f90d2407"],
Cell[CellGroupData[{
Cell[23536, 585, 631, 16, 72, "Subsection",ExpressionUUID->"e7d99775-53c8-443d-bc9b-999bb355f654",
 Evaluatable->False],
Cell[24170, 603, 1821, 50, 375, "Input",ExpressionUUID->"bbb28080-82de-4f03-b226-4247ac15aab6"]
}, Open  ]],
Cell[CellGroupData[{
Cell[26028, 658, 928, 20, 104, "Subsection",ExpressionUUID->"e57fe7c6-2316-45d0-9a5d-1efa01ac6413",
 Evaluatable->False],
Cell[26959, 680, 3459, 90, 375, "Input",ExpressionUUID->"157026c9-04c9-4d7a-9418-000cd75bd159"],
Cell[30421, 772, 148, 2, 47, "Input",ExpressionUUID->"dc939f53-588f-4af3-bd3c-17fb1dcbb632"]
}, Open  ]],
Cell[CellGroupData[{
Cell[30606, 779, 1224, 24, 104, "Subsection",ExpressionUUID->"391a5db5-3c8f-43f7-a11a-6166c9167d30",
 Evaluatable->False],
Cell[31833, 805, 744, 20, 90, "Text",ExpressionUUID->"8332d20d-e726-4921-9ed8-ebe81b7dc649"],
Cell[32580, 827, 3573, 92, 375, "Input",ExpressionUUID->"c68b342b-e3ef-4cb7-b2ab-44b43d21661b"],
Cell[36156, 921, 185, 3, 54, "Text",ExpressionUUID->"0b8915e2-a79f-4ff6-b058-57f0ca2c773f"],
Cell[36344, 926, 2486, 54, 345, "Input",ExpressionUUID->"26e627c5-dd72-4926-9caa-30c7417a47a4"],
Cell[38833, 982, 247, 6, 54, "Text",ExpressionUUID->"60d2d6fc-8561-4cd8-886b-50aad0651079"],
Cell[39083, 990, 409, 8, 54, "Text",ExpressionUUID->"607be239-0850-4601-b0da-2ef9ae900605"],
Cell[39495, 1000, 1732, 25, 44, "Input",ExpressionUUID->"ad49dce0-a7fe-4277-9a31-c57ad43d9c4e"],
Cell[41230, 1027, 1310, 28, 82, "Input",ExpressionUUID->"4dcd7886-c79d-4389-bfcd-3d40fd8c3417"],
Cell[42543, 1057, 449, 9, 47, "Input",ExpressionUUID->"c4fb0e91-6165-42ff-bd0a-16a892ae8a00"],
Cell[42995, 1068, 3115, 81, 524, "Input",ExpressionUUID->"8c2a8527-99b1-4999-8ec7-13a9dff68758"],
Cell[46113, 1151, 338, 8, 47, "Input",ExpressionUUID->"732e3690-6c47-4bc2-a3e3-399ec1db597b"],
Cell[46454, 1161, 592, 13, 47, "Input",ExpressionUUID->"86949b1a-a94f-40fe-a008-106fcfaeae4d"],
Cell[47049, 1176, 197, 3, 47, "Input",ExpressionUUID->"55d9fa88-0013-43f6-85eb-5f40a1c1baa0"],
Cell[47249, 1181, 195, 3, 47, "Input",ExpressionUUID->"3d954a82-d59d-4213-ad29-836fb4c4a2b0"],
Cell[47447, 1186, 837, 22, 47, "Input",ExpressionUUID->"c2d1a6fc-b0e3-49c3-b5af-9d3ed52b0f86"],
Cell[48287, 1210, 999, 26, 47, "Input",ExpressionUUID->"513f85a1-1933-4b28-a727-4c13f8e32bc1"],
Cell[49289, 1238, 233, 5, 47, "Input",ExpressionUUID->"aabe84ee-187b-4050-aaae-51d1d8188fad"],
Cell[49525, 1245, 685, 14, 47, "Input",ExpressionUUID->"1d3a1415-ca3a-4501-8722-26383d7543ed"],
Cell[50213, 1261, 473, 12, 47, "Input",ExpressionUUID->"c1fbf734-e9ce-4f40-8932-6aadfcdcc4de"],
Cell[50689, 1275, 536, 16, 47, "Input",ExpressionUUID->"f06ae692-8295-4c7e-905a-e4d5cea4a4aa"],
Cell[51228, 1293, 219, 5, 47, "Input",ExpressionUUID->"14472c62-f4dc-4610-b5bd-a934de995d72"],
Cell[51450, 1300, 201, 3, 47, "Input",ExpressionUUID->"9d7133bb-c1d5-4c68-b218-809bad881501"],
Cell[51654, 1305, 195, 3, 47, "Input",ExpressionUUID->"84e4be1e-158f-4191-acf5-e481203f979e"],
Cell[51852, 1310, 320, 6, 47, "Input",ExpressionUUID->"9f9d6631-b1cb-4326-bfa6-c7d732cf2f38"]
}, Open  ]],
Cell[CellGroupData[{
Cell[52209, 1321, 1646, 42, 135, "Subsection",ExpressionUUID->"d92c15d0-00fb-4ff4-aabb-3632c6ddb818",
 Evaluatable->False],
Cell[53858, 1365, 247, 6, 54, "Text",ExpressionUUID->"2868ac19-adf4-4edc-92df-e570bddf8523"],
Cell[54108, 1373, 533, 17, 54, "Text",ExpressionUUID->"e3b11540-24d8-450f-9c78-55724d87f143"],
Cell[54644, 1392, 1793, 27, 44, "Input",ExpressionUUID->"6731fb10-3492-4361-bae2-795baac49b00"],
Cell[56440, 1421, 1340, 30, 82, "Input",ExpressionUUID->"8ec0362c-72d0-430d-bb9f-26e15a26d490"],
Cell[57783, 1453, 500, 10, 47, "Input",ExpressionUUID->"1fc4292c-3ccd-4980-a1db-09e0e4ea8d77"],
Cell[58286, 1465, 338, 8, 47, "Input",ExpressionUUID->"63a0e785-9b67-497d-b91a-0bd02737f2ac"],
Cell[58627, 1475, 635, 13, 47, "Input",ExpressionUUID->"3ac44665-b330-411c-8301-435d36908828"],
Cell[59265, 1490, 197, 3, 47, "Input",ExpressionUUID->"f7cf8411-0145-4367-9805-85220d4a15e9"],
Cell[59465, 1495, 195, 3, 47, "Input",ExpressionUUID->"bce37c15-84c9-4591-93a8-54fd7445d104"],
Cell[59663, 1500, 242, 8, 54, "Text",ExpressionUUID->"2c9309f8-1597-4f85-8bef-5d9f88aacc04"],
Cell[59908, 1510, 777, 21, 47, "Input",ExpressionUUID->"0a116422-d0d9-41a9-bdc2-f938ef32d89a"],
Cell[60688, 1533, 1195, 30, 47, "Input",ExpressionUUID->"1617cc33-a1a6-4c80-bf98-d6a813d7eb14"],
Cell[61886, 1565, 188, 4, 47, "Input",ExpressionUUID->"59f083c4-7601-459c-b5b3-93e7072b5f2f"],
Cell[62077, 1571, 633, 13, 47, "Input",ExpressionUUID->"1c49f1f1-c13e-4999-9be6-e7336c8513c2"],
Cell[62713, 1586, 516, 12, 47, "Input",ExpressionUUID->"b6e9f689-1e7a-4b7a-9e57-79620e8bf89e"],
Cell[63232, 1600, 261, 4, 47, "Input",ExpressionUUID->"ed5fd72c-853d-424a-91b8-99f5d2a46e1c"],
Cell[63496, 1606, 218, 5, 54, "Text",ExpressionUUID->"fcded526-a732-47b0-a93d-3bc366341383"],
Cell[63717, 1613, 710, 21, 148, "Input",ExpressionUUID->"c629bbb1-11a0-488c-8833-6bb096cc68aa"],
Cell[64430, 1636, 218, 3, 54, "Text",ExpressionUUID->"664e13cf-25b1-448c-bf9f-83c8c916553b"],
Cell[64651, 1641, 1056, 28, 148, "Input",ExpressionUUID->"6001bdcc-5a9a-4b01-91a9-db19fb9171d8"],
Cell[65710, 1671, 1378, 26, 409, "Text",ExpressionUUID->"efc5cf6b-ff48-4f83-833c-24a863cac639"],
Cell[67091, 1699, 501, 13, 86, "Input",ExpressionUUID->"a270da95-89f6-4280-b600-e6a07cf58d65"],
Cell[67595, 1714, 454, 13, 86, "Input",ExpressionUUID->"07ada612-cb0d-4a82-9cbb-d9fb2b98a5af"],
Cell[68052, 1729, 2519, 76, 109, "Text",ExpressionUUID->"6935b142-1c1c-464b-9871-f12ac94de75d"],
Cell[70574, 1807, 236, 8, 54, "Text",ExpressionUUID->"b66c14d1-2065-4d85-acac-4a5dac7a8f6b"],
Cell[70813, 1817, 533, 16, 47, "Input",ExpressionUUID->"4b636c6e-465d-4fe8-8d45-ed4f40e8522c"],
Cell[71349, 1835, 176, 4, 47, "Input",ExpressionUUID->"1986ff60-91c1-4b18-b28f-508c409cc8c6"],
Cell[71528, 1841, 201, 3, 47, "Input",ExpressionUUID->"1a6ae233-3ef8-404f-acbf-78216b4b1db1"],
Cell[71732, 1846, 195, 3, 47, "Input",ExpressionUUID->"92f91791-ba56-469f-b7d6-5dd728ceb56e"],
Cell[71930, 1851, 354, 7, 47, "Input",ExpressionUUID->"bea3a17f-9df5-4f0b-9b42-6637438532df"]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}
]
*)

